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ABSTRACT 

We present a new survey of HCN (1-0) emission, a tracer of dense molecular gas, focused on the little- 
explored regime of normal star-forming galaxy disks. Combining HCN, CO, and infrared (IR) emission, we 
investigate the role of dense gas in star formation, finding systematic variations in both the apparent dense gas 
fraction (traced by the HCN-to-CO ratio) and the apparent star formation efficiency of dense gas (traced by the 
IR-to-HCN ratio). The latter may be unexpected, given the recent popularity of gas density threshold models 
to explain star formation scaling relations. Our survey used the IRAM 30-m telescope to observe HCN(l-O), 

CO (1-0), and several other emission lines across 29 nearby disk galaxies whose CO (2-1) emission has previ¬ 
ously been mapped by the HERACLES survey. We detected HCN in 48 out of 62 observed positions. Because 
our observations achieve a typical resolution of ~ 1.5 kpc and span a range of galaxies and galactocentric radii 
(56% lie at > 1 kpc), we are able to investigate the properties of the dense gas as a function of local 
conditions in a galaxy disk. We focus on how the ratios IR-to-CO, HCN-to-CO, and IR-to-HCN (observational 
cognates of the star formation efficiency, dense gas fraction, and dense gas star formation efficiency) depend 
on the stellar surface density, Sstar, and the molecular-to-atomic gas ratio, Stnoi/Satom- The HCN-to-CO ratio is 
low, often ~ 1 /30, and correlates tightly with both the molecular-to-atomic ratio and the stellar mass surface 
density across a range of 2.1 dex (factor of 125) in both parameters. Thus for the assumption of hxed CO-to- 
H 2 and HCN-to-dense gas conversion factors, the dense gas fraction depends strongly on location in the disk, 
being higher in the high surface density, highly molecular parts of galaxies. At the same time, the IR-to-HCN 
ratio (closely related to the star formation efficiency of dense molecular gas) decreases systematically with 
these same parameters and is ~ 6 - 8 times lower near galaxy centers than in the outer regions of the galaxy 
disks. For hxed conversion factors, these results are incompatible with a simple model in which star formation 
depends only on the amount of gas mass above some density threshold. Moreover, only a very specihc set of 
environment-dependent conversion factors can render our observations compatible with such a model. Whole 
cloud models, such as the theory of turbulence regulated star formation, do a better job of matching our obser¬ 
vations. We explore one such model in which variations in the Mach number driving many of the trends within 
galaxy disks while density contrasts drive the differences between disk and merging galaxies. 

Subject headings: galaxies: ISM - galaxies: star formation - ISM: molecules - ISM: structure - radio lines: 
galaxies 
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1. INTRODUCTION 

It is not yet fully understood how the properties of molecu¬ 
lar gas clouds affect their ability to form stars. A strong rela¬ 
tion between gas density^ (n) and star formation is expected, 
because denser gas pockets will be more prone to collapse 
and will collapse more quickly. A key method to test this 
expectation has been to assemble observations of molecular 
lines with different critical densities and to compare these ob¬ 
servations to tracers of recent star formation. Because such 
observations do not need to resolve individual clouds or star¬ 
forming cores, this approach may be the most practical way 
to systematically study the effect of gas density on star forma¬ 
tion over large parts of the universe. In this paper, we present 
new observations of the HCN (1-0) transition that cover large 
(~ kpc) parts of nearby galaxy disks and compare them to es¬ 
timates of the recent star formation rate and CO emission, a 
tracer for the molecular gas content (H 2 ). HCN(l-O), with an 
effective excitation density > 10'* - 10^ cm“^, has (along with 
HCO^ (1-0)) become the most common tracer of dense gas 
over large parts of galaxies. By contrasting it with the bulk of 
the molecular gas traced by CO emission and the rate of re¬ 
cent star formation, we aim to gain insight into how physical 
conditions internal to the molecular gas, especially density, 
affect the star formation process. 


Starting with the groundbreaking 

papers by Gao & 

Solomon (^ 

.004a|b]l, several studies (e.g. 

Gracia-Carpio et al. 


pnog 

Bussmann et al.[2008l[Juneau et al.|2009 [Garcfa- 

Buril 

lo et a 

2012|l have shown that the total star formation 


linearly) with the luminosity of molecular lines that have high 
excitation densities (e.g., HCN) than with bulk tracers of the 
molecular gas (e.g., CO). While low-7 CO emission tracks 
star formation near ly linearly in the di s ks of local star form- 
ing galaxies (e.g ., Bigiel et al.| 2008] Schruba et al. 2011 


and [Leroy et~ar |26l3b| hereafter |L1 3 | l, the relationship be 

comes nonlinear (multivalued or steeper than linear) when 
starburst regions and merging galaxies are considered (e.g., 
Genzel et al.||2010| [Daddi et al.||2010|l. 


Considering spirals 
and (ultra)luminous infrared galaxies (Ltir > 10* * L©; here¬ 


after (U)LIRGs) as a single population, the SFR is a s uper- 
linear function of th e bulk gas tracer CO. By contrast, [Gao 
|& Solomon] ( |2004b) l showed that total SFR scales almost lin- 
early with the HCN(l-O) luminosity across the entire range 
of galaxies. In a follow-up study, jWu et al.j (|2005|l 
eluded that the large parts of galaxies sampled by |G; 


jao & 


Solo mon] ( |2004b| ) and Galactic dense cores align along the 
HCN(l-0)--IR relation. On the assumption that line lu¬ 
minosities are proportional to the mass of emitting gas, these 
results imply that the SFR per unit total molecular mass is 
higher in the brighter (higher SFR) systems, while the SFR 
per unit dense molecular mass is nearly constant. 

Explanations of the distinct relationships between CO, re¬ 
cent SFR, and HCN can broadly be grouped into two classes. 
One group of explanations, which we refer to as “density 
threshold models”, broadly posit that the SFR is determined 
by the mass of molecular gas above a certain threshold den¬ 
sity, with gas above this density forming stars at an approxi- 
mately fixed rate everyw here (i.e., SFR/Mjense is ~ constant). 
[Gao & Solomon] (j2004b| l and jWu et al.| ( |2005Tl advocat e such 

a view, as do jLada et al.] ( l2010]]20121 l and ]Evans et al.1 ( |20141 l 

Throughout this paper “density” refers to volume density of H 2 
molecules unless otherwise specified. 


based on studies of clouds in the Solar Neighborhood. A com¬ 
peting class of models, which we refer to as “whole cloud” or 
“turbulence regulated models,” posit that the properties of the 
whole star-forming cloud affect the efficiency of star forma¬ 
tion. That is, the average density, Mach number, and other 
properties of the entire cloud set the local threshold for star 
formation, the fraction of gas above this threshold, and the 


speed with which star-forming gas collapses. Eederrath & 
Klessen (2012)1 synthesize seve ral such models; in this pa- 


per, we f ocus on c omp arisons to]Krumholz & McKee ( [2005] 
hereafter ]KMK05]l and]Krumholz & Thompson] (]200'7] here¬ 
after [KTO?]!, which considered such a “whole cloud” picture 
in the context of observations of molecular lines with varying 
critical density. In these models, both the SER per unit of total 
molecular mass and the fraction of dense gas increase with the 
average density within the clouds (h). The Mach (Ai) number 
also affects the density distribution and criteria for star forma¬ 
tion, with a high Ai broadening the probability distribution 
function of gas densities and increasing high density conden¬ 
sations, though its numeri cal effect varies somewhat with the 
details of the model used ( Eederrath & Klessen|20121 l. 

Characterizing the relation between dense gas and star for¬ 
mation is critical to distinguish between these models, and, 
in turn, to understand how larger scale conditions in galax¬ 
ies affect star formation. Eor example, several recent results 
have highlighted the potential importance of dynamical ef- 
fects in explaining observed variations in the SER per unit gas 
([Genzel et al.]]2010[ [Daddi et al.]]20T0j [Garcia-Burillo et ak] 


2012||Meidt et al.|2013[ ). Eor a full picture of how such large 
scale phenomena propagate into cloud-scale variations in the 
star formation process, we must first understand how physical 
conditions within the gas (which may be aff ected by large- 
scale conditions, e.g., see Hughes et al.|2013| l affect star for¬ 
mation in the cloud. 

Observations of a diverse sample of normal galaxy disks 
can help distinguish these models but such a data set has re¬ 
mained lacking. The pioneering surveys of dense gas in galax¬ 
ies mostly offer rough averages over the entire molecu lar/star¬ 
forming disks or values at bright galaxy centers (e.g., [Gao &| 


Solomon|2004a[ Gracia-Carpio et al.|2006 2008 Bussmann 


et al. 

2008 Juneau et al. 2009] 

Krips et al. 2008| Crocker et al. 

TUn 

Garcfa-Burillo et al.|20] 

The existing handful of ob- 


have been time-consuming and p rov ide only small-size sam- 
ples (e.g., Brouillet et al. 2005 and Buchbender et al.|[2013| 


in M 31 and M 33, respectively)! Observations comparing ex- 
tragalactic observations to individual clouds in the Milky Way 

( [\^ et al.|2005[[Heiderman et al.|201^[Lada et al.|2012) have 

proved enlightening, but the Galactic clouds studied sample a 
more restricted range of conditions than nearby galaxies and 
also suffer from small sample sizes and narrow perspective. 

To bridge these studies, we undertook a new HCN sur¬ 
vey that targeted large parts of a significant sample of galaxy 
disks. Erom 2008 until 2011, we used the IRAM 30-m tele¬ 
scope to observe HCN (1-0) emission from the disks of 29 
nearby star-forming galaxies. We drew our targets from the 
HERACLES survey, which constructed m aps o f CO(2-l) line 
emission from a set of 48 nearby galaxie s (LI3 first maps and 
survey presented by [Leroy et al.|2009 1 . With a typical reso¬ 
lution of ~ 1.5 kpc, each of our HCN observations blends 
together a large population of individual molecular clouds. 
We expect this averaging to attenuate the intrinsic variabil¬ 
ity and evolutionary effects v isible when studying individual 
star-forming complexes (e.g., Schruba et al.|2010|l. However 
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this resolution is still fine enough that we can isolate many 
local variables related to the physical state of the interstel¬ 
lar medium (ISM), for example stellar and gas surface den¬ 
sity, interstellar radiation field, and so on. A particular em¬ 
phasis of our survey is to study the disks (not just the cen¬ 
ters) of star-forming galaxies, and so our observations span 
a wide range of galactocentric radii, from galaxy centers out 
to ~ 75% of the optical radius (defined as the radius of the 
25 mag arcsec^^ isophote in the B-band). This allows us to 
probe a wide range of local conditions and, because HERA¬ 
CLES builds on many surveys at other wavelengths (see Sect. 
1^, the data exist to readily estimate these local conditions. 
In the end, our HCN survey provides context for cloud-based 
studies in the Milky Way and Local Group galaxies, while 
contrasting the starbursts, whole-galaxy integrals, and bright 
galaxy centers targeted in previous extragalactic studies (see 
also Bigiel et al, subm.). 

This paper is structured as follows. We describe our ob¬ 
servations and the ancillary data sets in Sect. In Sect. 
we put our data in context by comparing them to star for¬ 
mation scaling relations obtained from unresolved observa¬ 
tions of galaxies. In Sect.|^ we focus specifically on our data 
and investigate the systematic variations across galaxy disks 
of the star formation efficiencies and the dense gas fraction 
in the molecular clouds. In Sect.|^ we compare these trends 
with the predictions from density-threshold and turbulence- 
regulated models of star formation. In Sect. we study how 
the compatibility between our observations and these models 
depends on the assumed CO and HCN conversion factors. In 
Sect. 1^ we try to define a common scenario for resolved and 
unresmved observations. Sect. [^summarizes the paper. 

2. DATA AND PHYSICAL PARAMETERS 
2.1. Data 

2.1.1. A New Survey of HCN (1-0) Emission in Galaxy Disks 

Erom 2008 through 2011, we used the IRAM 30-m tele¬ 
scope to observe HCN(l-O) line emission at 62 positions in 
29 galaxies drawn from the HERACLES survey (Table 
Eor comparison to the overall molecular reservoir, we also 
observed the CO(l-O) line at 58 of those positions (recall 
that HERACLES observed CO (2-1)). All of our targets are 
star-forming galaxies (Ltir < 10*^ Lq, hereafter SE galax¬ 
ies); most of our targets are disk galaxies, though several are 
known starburst galaxies (e.g., we include M 82) and several 
of our central pointings cover nuclear starbursts. The point¬ 
ings were chosen to cover a wide range in environment prop¬ 
erties and star formation activity, while also maximizing the 
likelihood of detection. To this end, we targeted regions with 
relatively bright CO(2-l) emission in the HERACLES maps, 
picking a set of such regions that span a wide range in galac¬ 
tic radius (up to 75% of the optical radius; see maps in Ap¬ 
pendix]^. Our final set of HCN( 1-0) detections spans a range 
of 2.1 dex (factor of ^ 125) in both stellar surface density and 
molecular-to-atomic mass ratio, a range of 1.9 dex (factor of 

80) in molecular gas surface density, and a range of 1.7 dex 
(factor of ^ 50) in SLR surface density. 

The angular resolution of the IRAM 30-m at the frequency 
of HCN(l-O) is ~ 28", which sets our working resolution in 
this paper. The corresponding spatial resolution at the dis¬ 
tances of our targets ranges from 0.4 to 3.4 kpc (rightmost 
column in Table with an average value of « 1.5 kpc. 

Eor observations in January and December 2008 we used 
the old AB receivers at the 30-m telescope. In August 2009 


and June 2011 we used the EMIR receiver ( [Carter et al.|2012[ 
see Table |^. Whenever possible with regard to the extent 
of the emission and the observing schedule, we observed 
in wobbler-switching mode with a total ON-OEE throw of 
+240" in azimuth and a wobbling frequency of 0.5 Hz. The 
orientation of the wobbler throw direction is fixed in azimuth, 
so the throw in equatorial coordinates depends on the hour 
angle of the target. We used the HERACLES CO(2-l) maps 
to ensure that every position was observed only when the 
OLE positions were free of molecular line emission. Due to 
this constraint, wobbler-switching mode was not viable for 
a few positions because the throw would have intersected 
the galaxy. In those cases we adopted a position-switching 
scheme and selected a suitable OFE position with coordinates 
chosen based on the HERACLES maps. 

We checked the focus of the telescope on planets or bright 
quasars at the beginning of each session and then every few 
hours and, if relevant, at sunset and sunrise. Every ~ 1- 
1.5 hours, we corrected the telescope pointing using a point¬ 
like source close to the target galaxy. The magnitude of these 
corrections typically agreed with the nominal pointing accu¬ 
racy of the telescope (~ 2" rms). Every 15 minutes, we ob¬ 
tained a standard chopper wheel calibration, which we used 
to place the data on the antenna temperature scale (T*^. In 
order to convert to main beam temperature (Tmb), we adopted 
forward and beam efficiencies from the IRAM documenta- 
tior pjth at were up-to-date at the time of the observations (see 
Table[^. The expected uncertainties in the resulting gain cali¬ 
bration are < 10%. We checked this by observing well-known 
line-calibrators on a non-systematic basi s. Comparing these 
to the 30-m telescope spectra collected byjMauersberger et al. 


(|1989|l, we confirm this accuracy. As a crosscheck, we also 
compared our HCN(l-O) line intensities at fi ve galaxy cen¬ 
ters to those measured by [Krips et ak] ( |2008j l, also with the 
IRAM 30-m telescope. The measured intensities agree within 
< 10 %. 


We reduced the data using the CLASS package of the 
GILDAS software librar}p] The raw spectra taken with 
EMIR (since 2009 onwards) are ~4 GHz wide. For each line 
at each position, we fitted a baseline to a 0.5 GHz-wide part of 
the spectrum centered on the line frequency Doppler shifted 
by using the systemic velocity of the observed galaxy. For 
observations taken with the AB receivers (in 2008), we used 
the entire 0.5 GHz bandwidth when fitting the baseline. De¬ 
pending on the quality of the baseline, we fitted a polynomial 
baseline of degree one or three. In each case we avoided the 
velocity range of molecular emission from the galaxy known 
from the HERACLES maps. When necessary, we removed 
spurious spikes from the spectra. In the positions observed in 
2009 (~ 20% of the total), we also suppressed a systematic 
small-amplitude ripple in the HCN(l-O) spectra via flagging 
and linear interpolation in the Fourier transformed spectrum. 
For each line and position, we averaged all the spectra ex¬ 
cept those that showed egregious artifacts, poor baselines or 
abnormal noise levels to produce a final spectrum. 

From these reduced, averaged spectra, we calculate 
velocity-integrated intensities, expressed throughout the pa¬ 
per on the main beam temperature (Tmb) scale in units of 
K km s '. We derive these intensities summing over the ve¬ 
locity range known to contain molecular line emission from 
the CO (2-1) maps (the same region that we avoided in base- 


http://www.iram.es/IRAMES/mainWiki/Iram30mEfficiencies 

http://www.iram.fr/IRAMFR/GILDAS 
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Table 1 

Galaxies Observed in the HCN (1-0) Survey 


Name 

NGC 

^Aj2000 

{h : m \ s) 

77ecj2000 
: ") 

D 

(Mpc) 

Ris 

(') 

i 

n 

PA 

n 

res. 

(kpc) 

0628 

01:36:41.8 

+ 15:47:00 

7.2 

4.9 

1 

20 

0.98 

2146 

06:18:37.7 

+78:21:25 

12.8 

2.7 

54 

123 

1.74 

2403 

07:36:51.1 

+65:36:03 

3.2 

7.9 

63 

124 

0.43 

2798 

09:17:22.8 

+41:59:59 

24.7 

1.2 

85 

152 

3.35 

2903 

09:32:10.1 

+21:30:03 

8.9 

5.9 

65 

204 

1.21 

2976 

09:47:15.3 

+67:55:00 

3.6 

3.6 

65 

335 

0.49 

3034 

09:55:52.7 

+69:40:46 

3.9 

5.5 

77 

68 

0.53 

3049 

09:54:49.5 

+09:16:15 

19.2 

1.0 

58 

28 

2.61 

3077 

10:03:19.1 

+68:44:02 

3.8 

2.7 

46 

45 

0.52 

3184 

10:18:17.0 

+41:25:28 

11.8 

3.7 

16 

179 

1.60 

3198 

10:19:55.0 

+45:32:59 

14.1 

3.2 

72 

215 

1.91 

3351 

10:43:57.7 

+ 11:42:14 

9.3 

3.6 

41 

192 

1.26 

3521 

11:05:48.6 

-00:02:09 

11.2 

4.2 

73 

340 

1.52 

3627 

11:20:15.0 

+ 12:59:30 

9.4 

5.1 

62 

173 

1.28 

3938 

11:52:49.4 

+44:07:15 

17.9 

1.8 

14 

15 

2.43 

4254 

12:18:49.6 

+ 14:24:59 

14.4 

2.5 

32 

55 

1.95 

4321 

12:22:54.9 

+ 15:49:21 

14.3 

3.0 

30 

153 

1.94 

4536 

12:34:27.0 

+02:11:17 

14.5 

3.5 

59 

299 

1.97 

4569 

12:36:49.8 

+ 13:09:47 

9.9 

4.6 

66 

23 

1.34 

4579 

12:37:43.5 

+ 11:49:05 

16.4 

2.5 

39 

100 

2.23 

4631 

12:42:08.0 

+32:32:29 

9.0 

7.2 

86 

86 

1.22 

4725 

12:50:26.6 

+25:30:03 

11.9 

4.9 

54 

36 

1.62 

4736 

12:50:53.0 

+41:07:14 

4.7 

3.9 

41 

296 

0.64 

5055 

13:15:49.2 

+42:01:45 

7.9 

5.9 

59 

102 

1.07 

5194 

13:29:52.7 

+47:11:43 

7.9 

3.9 

20 

172 

1.07 

5457 

14:03:12.6 

+54:20:57 

6.7 

12.0 

18 

39 

0.91 

5713 

14:40:11.5 

-00:17:20 

21.4 

1.2 

48 

11 

2.91 

6946 

20:34:52.2 

+60:09:14 

6.8 

5.7 

33 

243 

0.92 

7331 

22:37:04.1 

+34:24:57 

14.5 

4.6 

76 

168 

1.97 


Note. — Center coordinates, distan ce, size, and orientation of the observed 
sources. All values are taken from |L12| when av ailable. The remainde r objects are 
NGC 2146 and NG C 2798 (all data take n from |Schruba et a l |201H, NGC 3034 
(all from HyperL eda.IPaturel et al.|2003| except distance, fr omlWalter et al. 2002^ , 
NGC 3077 (from |WaIteret^.|2uJs( , and NGC 4631 (from llrwmeTaT'PP' ^7 ' 


line fitting). We estimate the statistical error, crune, via 

tT'line “ ^ITiychannel V^VA V, 


( 1 ) 


where rmichannei is the rms noise, in Kelvin, for each channel, 
6V is the channel width and AV is the width of the line win¬ 
dow. Throughout this paper, we disregard the fine structure of 
HCN, which is not apparent in our observations, and quote a 
single integrated intensity for the 7 = 1 —> 0 transition. The 
two brightest fine-structure levels of the 7 = 1 —» 0 transition 
lie only ~ 5 km s ' apart. This is comparable to our typical 
spectral resolution (3.4 - 6.8 km s ') and three times smaller 
than the smallest line width (FWHM) of any detected line in 
our survey. 

We acquired and reduced the CO (1-0) data in a similar 
way, although the integration times were much shorter than 
for HCN (1-0) because CO (1-0) is a much brighter line. In 
most of this paper we derive molecular surfa ce densities from 
the CO(2-l) HERACLES maps (Sect. [ZS, because, unlike 
the CO(l-O) data , they are available at all observed positions. 
In Appendix [at] we use the CO(l-O) data where available to 
demonstrate that the CO excitation conditions do not bias our 
main results. 

Appe ndix |B] s ummarizes the results of our 30-m survey. 
In Figs. B1 -[B5| we show our 62 observed positions on the 


HER ACLES CO(2-l) maps of the target galaxies. Figs. B6 - 
|B19| sh ow the final HCN(l-O) and CO(l-O) spectra. Ta¬ 
ble B1 reports the velocity integrated intensities, /hcnio and 
^coio^ at each position, indicating statistical errors, cthcnio 
and cTcoio, in brackets. The quoted errors do not include the 

10% uncertainty that we expect in the overall flux calibra¬ 
tion. For the analysis in this paper, we consider all spectra 
with SNR\ine = (/line/cTiine) > 4 to be detections. As a check, 
we verified that in all cases the line shapes of HCN were con¬ 
sistent with the shape of the CO (2-l) line at matched spa¬ 
tial and spectral resolution (Figs. |B6f|B19|l. We applied the 
same check to spectra with lower signal-to-noise ratio and, 
based on this criterion, we promoted three HCN spectra with 
3.4 < SNRhcn ^ 3.7 to be detections. In summary, we detect 
HCN(l-O) at 48 positions with a median 5 ARhcn of ~ 9.5. 
The HCN detections span ~ 2.5 dex in /hcnio (a factor of 
^ 330) and extend from galaxy centers up to r 25 = 0.75. The 
detection at the largest radius occurs in NGC 6946 at a galac- 
tocentric radius of 8.5 kpc. We detect the CO(l-O) line in 57 
out of 58 positions with a median SNRco of ~ 35. When 
reporting non-detections, we give 4cr upper limits. 

2.1.2. Ancillary Data 

Our target galaxies belong to the HERACLES sample, 
which means that we have CO(2-l) emission mapped at 
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Table 2 

Receiver-Backend Combinations. 


Month/year 

Receiver(s) 

Backend 

01/2008 

AlOO-i-BlOO 

1 MHz filterbank 

12/2008 

AlOO-i-BlOO 

1 MHz filterbank 

08/2009 

E090 

WILMA (2 MHz resolution) 

06/2011 

E090 

WILMA (2 MHz resolution) 


Table 3 

Observed Lines 


Line 

Frequency 

Beam 

Tmb/T^ 


(GHz) 

(") 

[2008] [2009/2011] 

HCN(l-O) 

88.6 

27.8 

1.23 1.17 

CO(l-O) 

115.3 

21.3 

1.28 1.22 


Note. — Line rest frequencies. Reported beam size corre¬ 
sponds to the FWHM of the 30-m beam. Separate Lmb/T^ 
ratios are given for the 2008 and 2009/2011 campaigns, which 
used different receivers (see text and Tablej^for details). 


~ 13" across their disks. HERACLES targets many individu¬ 
ally well-studied nearby galaxies and was designed to overlap 
with other surveys that provide data across the electromag¬ 
netic spectrum. The data from these other surveys and from 
numerous studies of individual galaxies, provide us with an 
excellent characterization of the local properties of the ISM 
and of the star formation rates at each of our poin tings. A 
detailed account can be found in|Leroy et al.| ([20121 hereafter 
LI2; see also |L13||. In short, we make use of the following 
data; 


CO(2-l) line emission, from HERACLES ( |L13| ). 

IR emission in the 3.6 - 160 /tm range, from the 


Spitzer Infrared Galaxies Survey (SINGS, Kennicutt 
et al. 2003 i and t he Local Volume Legacy Survey (LVL, 
|Dale et ST 2009 [ ). 


Continuum-subtr acted Ha emission, from 

LVL, GoldMine ( |Gavazzi et al.||2003]l, and __ 

lation of other, smaller studies ([Hoopes et al. 2001 


SINGS 
a compi- 


Boselli & Gavazzi|2002[ Knapen et al.|2004]r 


Hi line em ission, from The H i Nearby Galaxy Survey 
(THINGS, [Walter et al.|200^ and a collection of new 
and archival VLA data (|L13 and Leroy et al., in prep.). 


Except in the few cases noted below, we convolved all the 
available maps to our working resolution of 2 8", i.e ., the reso¬ 
lution of our HCN data. We closely followed |L12| to estimate 
physical parameters from them. 


2 . 2 . Physical Parameter Estimates 

All surface densities in this paper art face-on values. Eor 
this, the equations to convert intensities measured along the 
line of sight to surface densities include a cos(0 factor, where 
i is the adopted inclination of the galaxy disks (Table[^. The 
cosine factor corrects for the observational bias that makes 
measured intensities to increase with i as the beam intersects 
a larger volume of the disk. 


In all cases we have tried to keep the physical parameters in 
this paper closely linked to observed intensities. That is, total 
molecular surface density is a recasting of CO intensity, dense 
gas surface density is a linear translation of HCN intensity, 
and so on for Hi (21-cm intensity), star formation (IR intensity 
from 24/im), and stellar surface density (near-IR intensity). 
We explore additional effects related to parameter estimation 
in the appendix but the results of this paper can mostly be 
straightforwardly read in terms of observables and we plot 
these as an alternative axis whenever possible. 


2 . 2 . 1 . Molecular Surface Density 

The mass surface density of molecular gas (Snjoi, including 
helium) is commonly derived from the velocity-integrated in¬ 
tensity of the CO (1-0) line. Eor practical reasons, in this pa¬ 
per we derive Emoi from the HERACLES observations of the 
CO (2-1) line, /co 2 i; we can easily convolve our full maps 
at higher resolution to match the HCN(l-O) beam and derive 
a CO(2-l) intensity for each point (while the CO 1-0 line is 
observed at higher resolution than the HCN 1-0 line because 
of its rest frequency). We then define Iqq as the CO(l-O) 
intensity inferred from the HERACLES observations on the 
assumption of a fixed CO(2-l)-to-CO(l-0) line ratio, RaL 


^co - 


JC021 

^21 


( 2 ) 


We a dopt R 21 = 0.7, which is the average ratio found 
by |L13|by comparing their HERACLES CO(2-l) data with 
CC)(l-0) data available in the literature (including the data 
presented in this study) at matched angular resolution. Our 
new observations of CO(l-O) support this char acteristic ratio 
as appropriate for our data. In Appendix |A.1| we redo large 
portions of our analysis with the available CO (1-0) data and 
confirm that the assumption of a fixed R 21 does not bias our 
results. 

Erom I CO, we derive 


^mol “ Otcolco cos(0, 


(3) 


where aco is the CO-to-molecular mass conversion factor, 
here phrased in terms of CO (1-0). By default, we as¬ 
sume throughout our sample a fixed. Galaxy-like conver¬ 
sion factor aco - Q'°q = 4.4 Mq pc“^ (K km s“')“' (de¬ 
rived from a H 2 column density conversion of Aco - 2 x 
10^*^ cm“^ (Kkm s“')“' and a 1.36 helium-correction factor). 
We explore the impact of varying this assumption in Sect.|^ 
Note that our approach to aco could simply be recast as 
adopting a fixed conversion factor for the CO (2-1) intensity. 
Eor an in-depth investigation into the CO-to-H 2 con version 
factor in HERACLES galaxies, we refer the reader to jSand- 
jstrom et ar] ( |2013| |. 


2 . 2 . 2 . Dense Molecular Gas Surface Density 

Molecular clouds host cores and dense filaments comprised 
of gas with densities several orders of magnitude higher than 
the cloud average (n). The emission lines of molecules with 
high electric dipoles (e.g., HCN, HCO^, CS) trace this dense 
component, in which collisions can thermally excite their en¬ 
ergy levels. A conversion factor (cthcn) can be defined to cal¬ 
culate the mass surface density of dense molecular gas (in¬ 
cluding helium), Sdense, from the HCN(l-O) intensity: 


^dense “ ancNfucNio cos(i). 


(4) 
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We initially 
10 Mq (K km 


assume a 


fixed 


ffHCN - 


[qi. 


Gao & Solomon (2004b I 


o (K Km which 

argued to be reasonable for the disks of normal, star-forming 
galaxies. In discussing Sdense, we consider this quantity to 
represent the surface density of gas at densities above a fixed 
cutoff, Udense, where Wdense IS of the order of the star formation 
threshold expected from some models and observations, i.e., 
~ 10“*-10^ cm“^. The exact value of Wdense does not matter for 
most of our analysis, as long as this density is associated with 
the dense, often star-forming substructures within clouds. 

A potentially important issue is how Wdense compares with 
the density required to effectively excite the HCN(l-O) line 
(Heft)- For the HCN(l-O) line, the critical density at which the 
de-excitation rates by collisions and by spontaneous emission 
are equal is Merit = 1.25 x 10® cm~^ at Tk = 20 K (taking 


the rate coefficients listed in the LAMDA databasd'^ Schoier 


|et al.|2005] l. However, it is well-known f rom radiative transfer 
models (e.g.,|Scoville & Solomon] 1974|l that a line can be ef¬ 
fectively excited at densities lower than Merit thanks to opacity 
effects. The excitation threshold is thus sensitive to variations 
in the gas temperature, the velocity gradient, and the HCN 
abundance within the clouds. Based on a grid of models with 
the RADEX radiative transfer code ( |van der Tak et al.|2007| l, 
we hnd n^ff could be ~ 10^ - 10® cm“^ for typical cloud con- 
ditions. With our dehnitions, any variation in this excitation 
threshold could imply a change in othcn- We explore how an 
Q'hcn that varies with environment would impact our conclu¬ 
sions in Sect.|^ 

2.2.3. Other Local Conditions Within Our Beams 

Atomic gas. We assume that the 21 cm Hi emission is op¬ 
tically thin and derive the mass surface density of atomic gas 
(Zatom, including helium), from the velocity-integrated inten¬ 
sity of the line, /hi: 


Mopc- 


= 1.98 X 10" 




K km s 


■ cos(/). 


(5) 


The total mass surface density of gas, Zgas, is defined as 
2^gas = 2inioi + ^^atom and the molecular-to-atomic ratio is 
2^moi/5^atom- In the observed positions where HCN is detected, 
the molecular gas makes up ~ 50 - 100% (median 88%) of 
the total gas surface density. 

Stars. We estimate the mass surface density of old stars 
from the 3.6 fim intensity, /3 6 via; 


Mo pc 


= 280- 


h .6 


MJy sr 


q- cos(0, 


(6) 


which we consider uncertain by ~ 50% due to contamina¬ 


tion from young stars and dust emission ( Meidt et al. 2012 
and uncertainties in the mass-to-light ratio (|Zibetti et al. 2009 
IMeidt et al.|M4l ). 

Dust Properties and IR Spectral Energy Distribution. We 
convolve the LVL and SINGS Spitzer images to a common 
resolution (~ 40"), construct radial prohles for each galaxy, 
and fit the average spect ral energy dis t ributio n in each ring 
using the dust models of [Draine & Li|(|2007|l; for details of 
this specific processing see [Leroy et al. (2012 1 . From these 
models we derive several dust properties within radial bins, 
of which the 24 /rm-to-TIR ratio, the “cirrus” emissivity, and 
the dust-to-gas mass ratio (after comparing with CO and Hi 

** http://home.strw.leidenuniv.nl/~moldata/ 


data also convolved to 40") are used in this p aper t o derive 
star formation rates (see below and Appendix |A.2|i. These 
quantities are all ratios (rather than absolute intensities) and 
we assume that they hold, at least on average, at the somewhat 
higher working resolution of our HCN observations. 

Star Formation Rate (SFR). To easily compare our results 
with most studies of molecular lines in galaxies, we estimate 
the SFR surface density (Zsfr) from the total IR intensity 
(Itir); 


^SFR 


Mq Myr^' pc“ 


= 1.87X 10"® 


friR 


Lq pc 2 sr 


cos(0. (7) 


Here, /tir is calculated from the measured 24 pm in¬ 
tensity (I 24 ) and the dust models, which produce a 24-to- 
TIR ratio that rigorously holds at slightly larger scales than 
our beam. This can alternatively be thought of as approxi¬ 
mately using the 2Apm emission to make an aperture correc¬ 
tion that matches the lower resolution TIR, estimated from 
24 pm, 70 pm, and 160 pm data, to our 28" H CN beam. 

Equ ation [7| c orresponds to the calibratio n by [Murphy et al.j 
( [201 1[ see ^o [Kennicutt & Evans||20T^ and we prefer this 
approach because it links us cleanly to previous work compar¬ 
ing HCN and CO, which has focused mostly on the integrated 
IR luminosity as a star formation rate tracer. This allows our 
results to be seamlessly read in terms of IR luminosity su rface 
density, even though we label them “SFR.” In Appendix [A.2[ 
we take advantage of the ancillary data set and build alterna¬ 
tive SFR tracers (Ha, 24/im, and a linear combination of both) 
to confirm that the choice of star formation rate tracer has a 
minimal impact on our results. A ll the SFR calibrations in 
this paper assume a Kroupa (2001[l IMF. 


2.3. Extragalactic Reference Sample 

We compare our data to observations of mo lecular lines in 
106 SF galaxies and (U)LIRGs compiled by |Garcla-Burillo[ 
et al. ( 2012[ hereafter GB12). GB12 combined published 


data with new observations using the IRAM 30-m telescope. 
GB12 present CO(l-O) and HCN(l-O) luminosities for all 
of their targets, estimating these from higher-order transi¬ 
tions in a few high-redshift objects. Their sample includes 
both upper limits (non-detections) and lower limits for ob¬ 
jects without full mapping. In principle, we derive the total 
SFR (in Mq Myr ') for these targets from the total IR lu- 
minosity (Ltir) again using Equation |7] However, following 


Gracia-Carpio et al. ( [2008^ , we rather use 1.3 x Lfir instead 
of Ltir. where Lfir is the far IR luminosity. This allows us 
to avoid AGN-driven contamination of Ltir in the brightest 
(U)LIRGs, while it hardly changes the estimated SFR in ob¬ 
jects powered by star formation. In particular, once the most 
conspicuous AGN sources are removed, the Ltir/Tfir ratios 
of the GB12 sample show a small scatter of 0.05 dex rms 
(12%) around a median value of 1.3. Since we lack TIR 
data for some GB12 objects, we also use (1.3 x Lfir) rather 
than Ltir to classify them as either SF galaxies or (U)LIRGs 
(Sect. [T]|. In principle, this affects the classiheation of only 
hve galaxies. 

For simplicity and in order to present the data as directly as 
possibl e, w e assume the constant conversion factors quoted 
in Sect. [2.2[ to derive (dense) molecular masses from the mea¬ 
sured line luminosities. This is a non-trivial assumption, given 
the ample evidence that these factors are a few times lower in 
(U)LIRGs than in SF galaxies (e.g., see Bolatto et al.||201? 
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for a review). The main influence of these factors on our first- 
order analysis is likely to be some offset between (U)LlRGs 
and SF galaxies. We return to the topic of conversion fac¬ 
tors in Sect. An additional concern is the unobscured SFR, 
which is not perfectly accounted for by the IR luminosity and 
could make up a significant fraction of the tot al SFR in SF 
galaxies, thus causing additional offsets (|GB12)l. 

For the 83 objects that have measured sizes for the star¬ 
forming/molecular disk, we convert luminosity/flux measure¬ 
ments to estimates of the average £den.se, ^inioh and Ssfr- We 
divide the total molecular (dense) mass or SFR in each galaxy 
by its disk area. The measured sizes are typically given 
as FWHM values of the light distribution, so that we quote 
(2) = M/(7rFWHM^^j/4). Here, FWHM„,aj would be the 
Full Width at Half Maximum along the the major axis of the 
star-forming/molecular disk. This definition is equivalent to 
that o f the face-on surface densities derived from our data 
(Sect. |2.2| i. When calculating (S), we implicitly neglect any 
aperture correction to account for the fraction of the mass in¬ 
side the FWHM, which would be a factor of 0.5 for a two- 
dimensional Gaussian and potentially as low as ~ 0.15 for an 
exponential disk. It is also driven by the total mass and so 
may be viewed as somewhat closer to the peak S value than 
the value likely to be found at any given locale in the galaxy 
(in the same way that the ratio between the total area below a 
Gaussian curve and its FWHM approximately gives its peak 
value). For the example of a face-on exponential disk, the 
surface density (S) is, in fact, ~ 4 times the peak value of 
the exponential because the aperture correction (i.e., the mass 
outside the FWHM) is so large. For a Gaussian, the number 
approaches ~ 1 .4. Most SF galax ies likely fall somewhere in 
this range (e.g.,|Young et al.|1995|l. We return to this source of 
bias when we compare our sampling of galaxy disks to these 
unresolved data (Sect. [^. For now, we caution the reader that 
(L) derived in this way is biased somewhat high. 


2.4. Analysis tools 

We parameterize the relations between pairs of variables 
using power-law fits, which we use to capture general trends 
in the data. We derive the best fit power law from total-least- 
square minimization, typically in the log-log space, excluding 
non-detections and assuming equal uncertainties along both 
axes. Quoted uncertainties in the fit outputs are Icr levels ana¬ 
lytically derived from the equations of error propagation. Our 
lower/upper limits are typically compatible with the solutions 
that we find, so we do not expect that their exclusion biases 
the fit outputs. Given the irregular sampling of our targets, we 
caution that care should be taken before extrapolating our fits 
to other regimes or data sets. 

We also find it useful to measure the strength of the corre¬ 
lations between pairs of variables in a non-parametric way, 
which we do using Spearman’s rank correlation coefficient 
(rs; hereafter rank coefficient). Throughout this paper, rank 
coefficients do not consider limits (i.e., non-detections). We 
measure the significance of by its p-value, i.e., the proba¬ 
bility that the absolute value of the rank coefficient for uncor¬ 
related data could be equal or higher than the measured value 
by chance. Throughout this paper, we consider that a corre¬ 
lation is significant when the p-value of its rank coefficient is 
lower than 2.5%. 

We performed most of the data analysis within the environ- 
ment of statistical computing R (|R Development Core Team 
[20T3] l. 


Table 4 

Rank Correlations for (Dense) Molecular — SFR Scalings 



Mdense 

SFR 

Mmol vs. 
SFR 

^dense VS. 

2^sfr 

^mol VS. 
^SFR 

this paper+QE\2 

0.96* 

0.96* 

0.97* 

0.97* 

this paper 

0.89* 

0.95* 

0.88* 

0.87* 

GB12 

0.89* 

0.86* 

0.92* 

0.89* 

GB12-SFonly 

0.54* 

0.55* 

0.95* 

0.92* 

GB12-(U)LIRGonly 

0.82* 

0.83* 

0.91* 

0.89* 


Note. — Rank correlation coefficients relating star formation and gas sur¬ 
face densities (Fig.[TJ for various groupings of our data and the GB12 data. 
Star symbols indicate significant correlations (p-value < 2.5%). 


3. COMPARISON WITH UNRESOLVED OBSERVATIONS: STAR 
FORMATION SCALING RELATIONS 


Our HCN detections span a range in HCN (1-0) luminos¬ 
ity (within the beam) of ~ 3 x 10^ - 2 x 10^ K km s“^ pc^. 
Thus, they span the luminosity gap betwe en individual Ga lac¬ 


1 xlO"^ Kkm s ' pc^ in 


Wu et al. 


T- 


2005 I and 


GB12I. 


tic clouds (up to 

entire galaxy disks (down to ~ 1 x 10' Kkm s ‘ pc" in 
This reflects that, at our working resolution of ~ 1.5 Kpc, we 
probe ensembles of clouds at sub-galactic scales. Because this 
scale has been barely explored in previous studies, we begin 
our analysis by comparing in Fig. [T] our data (gray circles) to 
unresolved observation s of galaxies drawn from GB12 (col¬ 
ored symbols; see Sec t.|2.3| l. 

The top row of Fig. [H plots the total SFR in our pointings 
(within the 28" beam) and in the GB12 galaxies as a func¬ 
tion of the (dense) molecular mass traced by HCN (Mdensei 
left-hand column) and CO (Mmoi; right-hand column). The 
bottom row shows the same information, but in terms of sur¬ 
face densities of SFR and (dense) gas mass. These surface 
density scaling relations are often called Kennicutt-Schmidt 
(KS) plots, and the accompanying power-law fits (“KS laws”) 
are a common tool for empirical studies of the relations be¬ 
tween star formation and the ISM. We remind that, for sake 
of simplicity, we adopt fixed CO and HCN conversion fac¬ 
tors for both SF galaxies and (U)LIRG. Thus, these are also 
luminosity-luminosity scaling relations and we include axes 
showing the corresponding IR and line luminosities. 

As expected from previous studies, in all four panels of 
Fig .[^the SFR shows a good correspondence with the (dense) 
molecular gas traced by the HCN and CO emission. We also 
see this in Table which lists the rank coefficients for the 
plots in Fig. considering various groupings of our data and 
the GB12 data. All of these correlations are significant (p- 
value < 2.5%) and high (> 0.8 for total SFR/masses, except 
for the subset of SF galaxies from GB12; > 0.9 for surface 
densities), thus reflecting the tight correlations seen in the 
figures. Unexpectedly, the rank coefficients of the SFR-vs- 
^dense SFR-vs-M^oi relations, on the one hand, and of 
the SsFR-vs-Smoi and SspR-vs-Sdense relations, on the other, are 
very similar for every data set. Although the slopes of the re¬ 
lations vary, implying differences in the physics at play for 
different molecular tracers, the CO and HCN lines appears to 
be almost equally good at predicting the SFR at these scales. 
We caution the reader that this result may not hold at different 
scales and/or different data samples. 

Fig. n shows that our data extend the range of gas masses 
and surface densities spanned by unresolved observations. 
In terms of total gas masses (top row), our brightest points. 
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Figure 1. Top row: Total SFR (Ssfr) as a function of the mass of molecular (left) and dense (right) molecular gas for pointings in galaxy disks (this paper, 
gray points) and unresolved star-forming galaxies (red) and (U)LIRGs (blue) from GB12. Bottom row: Surface density of recent SFR (I^sfr) as a function of 
the surface density of total (left) and dense (right) molecular gas for the same data sets. The top and right-hand axes of each panel display the data in terms of 
observed quantities. Open symbols indicate limits in the direction of the attached arrows. The equations in the figures report power-law fits to our data (gray 
line), the GB12 sample (green line), and all the data (brown line). Eri'ors at Icr-level in the fit parameters are indicated. For comparison, the yellow area shows a 
fixed ratio (power law index 1) with a factor of 2 scatter. 


mostly vigorously star-forming galaxy centers, lie right below 
the least active SF galaxies (dark red circles). In terms of sur¬ 
face densities (bottom row), there is a partial overlap between 
them. Pointings at larger galactocentric radii, where the lines 
are typically fainter (Sect. |^, have masses and surface den¬ 
sities as much as two and one order of magnitude lower than 
the galaxies in the GB12 sample, respectively. Recall that 
we expect the surface densities estimated from the unresolved 
(GB12) observations to be bias ed s omewhat high by the lack 
of an aperture correction (Sect. [23]l. Correcting by these fac¬ 


tors, which should affect both axes, would tend to increase the 
GB12 overlap with our data set in terms of surface density. 

The bottom-row panels of Fig.[2might support the idea that 
(U)LlRG s and SF galaxies form distinct star formati on se¬ 
quences ( |Genzel et al.|2010]|Daddi et al.|2010]|GB12| l. With 
or without aperture correction applied to the unresolved ob¬ 
servations, our data and the SF galaxies from GB12 align 
along a continuous and monotonous sequence. In contrast, the 
(U)LIRG sample lies parallel, but seemingly offset to higher 
2isFR for the same (dense) gas surface density. The offset is 
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particularly noticeable in the Ssfr - Smoi plane, and would be 
even larger if we assumed lower-than-Galactic CO and HCN 


conversion factors for the (U)LIRG subsample of GB12 

In the four panels of Fig. adding our data to the GB12 
sample has the effect of lowering the best-fit power law in¬ 
dices. However, the reason for this change differs between 
the left-hand and right-hand panels. In the Ssfr - Smoi plane 
(similarly in the SFR - M^o\ plane), the global power-law 
index drops from ~ 1.4 to ~ 1.2 because our data (power- 
law index of 1.01 + 0.08) reinforce the linear trend shown by 
the SF galaxies in the GB12 sample (index of 0.97 + 0.07; fit 
not shown in the plots). By contrast, in the Ssfr “ ^ijense plane 
(similarly in the SFR- plane) the global index changes 
from ~ 1.1 to ~ 1.0 when our data are included because our 
data by themselves follow a significantly sublinear relation 
(index of 0.69 ± 0.06). In this plot, the SF galaxies alone are 
fitted by a higher index of 0.88 + 0.06. This highlights a main 
conclusion of this paper; that, perhaps surprisingly, the ap¬ 
parent efficiency with which dense gas forms stars seems to 
vary systematically across our sample of galaxy disk point¬ 
ings. The sense of the variation is that regions with lower 
dense gas fractions and lower dense gas surface densities tend 
to have a higher apparent rate of star formation per unit dense 
gas. 

In summary, despite the continuity between our data and 
SF galaxies. Fig. [T] suggests that the star formation laws for 
the dense molecuK gas are not the same at galactic and sub- 
galactic scales. This is not illogical, given that many of the 
new conclusions in this paper come from data points with low 
surface density. These typically correspond to off-center po¬ 
sitions and, because they are relatively faint in line emission, 
they will not dominate the total galaxy luminosity. Our re¬ 
sults thus suggest that the average star formation properties of 
molecular clouds vary across galaxy disks. We investigate this 
topic in detail in the following sections, focusing on our new 
data set. Doing so, we take advantage of the fact that our data 
set resolves disk galaxies into discrete regions and targets the 
well-characterized HERACLES (i.e., THINGS and SINGS) 
galaxies. 

4. DENSE GAS ERACTION AND STAR FORMATION EFEICIENCY 
ACROSS GALAXY DISKS 

4.1. Normalized Quantities: /denser SFE^^i, and S FE^^nse 

In this section, we combine our HCN obser vation s with 
ancillary data of the observed pointings (Sect. |2.1.2| i. We 
explore how the beam-averaged properties of the dense gas 
in unresolved molecular clouds depend on local conditions 
within a galaxy. Doing so, it is crucial to bear in mind that, 
because the molecular interstellar medium in galaxy disks 
can be sparse and clumpy at our typical spatial resolution, 
neither beam-integrated masses (“luminosities”) nor beam- 
averaged surface densities (“intensities”) can be straightfor¬ 
wardly transl ated into intrinsic clo ud properties (see extended 
discussion in Leroy et^|2013^ . For example, in the sim¬ 
plistic assumption that the beam had an area A and probed N 
identical spherical clouds of mass mdoud and projected area 
Qcioud, we would find that 


— nicioud ^ ^ ? 


(8) 


N 

^ — ^cloud ^ ^cloud ^ 


(9) 


where Sdoud = iWcioud/flcioud is the intrinsic surface density in a 


cloud and M and Z are the mass and surface density measured 
at beam scaleJ^ 

Because we work at a fixed angular resolution, beam- 
integrated masses (M) are strongly dependent on the source 
distance, that determines the number of clouds captured by 
the beam {N) to a significant extent. This observational 
bias is partially corrected when working with surface densi¬ 
ties. However, because our observations span a large range 
of galactocentric radii, we expect large variations in the NjA 
ratio, which measures the “filling factor” of clouds per unit 
area. These filling factor variations, which may largely affect 
SsFR, S^oi, and Sdense in the same way, can drive surface den¬ 
sity correlations (e.g., see Sect.j^and Table|^, masking more 
subtle changes in intrinsic cloudproperties. 

The simplest way to eliminate the Af/A factor, assuming 
that it acts approximately equally on all surface density terms, 
is to represent our data in terms of three surface density ratios 
that are also frequently used in the literature; the dense gas 
fraction (/dense = 2dense/2imoi), the Star formation efficiency of 
the molecular gas {SFE^o\ = SsFR/^moh i-S-, the inverse of 
the molecular gas depletion time, with dimension of time *), 
and the star formation efficiency of the dense molecular gas 
(SFE^gnse = ^SFR/2idense)- These three parameters are related 
through the following equation; 


^ FE[hqI — S EEdense ^ /dense ■ (10) 

If all the star formation and molecular emission arises from 
the same set of clouds inside a beam, then these ratios cap¬ 
ture the average intrinsic cloud properties over the population 
within the beam, independent of any filling factor variations. 

Our simple approach to physical parameter estimation 
means that for most of this paper these rations have direct ob¬ 
servational analogs. The ratio of HCN-to-CO emission maps 
to the dense gas fraction, the ratio of infrared to HCN emis¬ 
sion maps directly to the dense gas star formation efficiency, 
and the ratio of infrared to CO emission maps to the overall 
molecular gas star formation efficiency. 

4.2. Environmental Trends 

After exploring a variety of ISM parameters, we found 
that the systematic variations of /dense, SFE^ense, and SFE^oj 
in our sample are best captured by correlations with the 
molecular-to-atomic ratio, Smoi/Satom, and the stellar mass 
surface density, Sstar- These two quantities also have the ad¬ 
vantage of being independent of one another and representing 
extremely simple functions of the available data, which allows 
us to confirm with little analysis that those correlations are not 
spurious. We stress, however, that many other ISM parame¬ 
ters that are relevant to our study are highly covariant with 
2imoi/2atom and Sstar (c-g-, metallicity, midplane pressure), so 
that we cannot rigorously ascertain what ultimately drives the 
trends discussed below. The reader should keep in mind that 
both Smoi/diatom and Sstar decline with increasing galactocen¬ 
tric radius, so that the lower values in Fig. [^typically corre¬ 
spond to outer positions in the disks. 

Fig. 1^ shows /dense IhCn/IcO'- tOp rOW), ^EEdense (°^ 
Ar/^hcn; middle row) and SFEj^oi ^ir/^co^ bottom row) 
as a function of Smoi/^atom (left column) and Sstar (right col¬ 
umn). The Figure shows several well-defined, systematic 

As usual in extragalactic studies, we assume that the cloud-cloud veloc¬ 
ity dispersion prevents clouds in the foreground from absorbing the emission 
of those behind them. 
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Table 5 

Rank coefficients for /denser S and 

S -FEfnoi as a function of local conditions in the 
galaxy 



/dense 

S FEiensa 

S Ffimoi 

'"25 

-0.59* 

0.43* 

-0.07 

^mol/^atom 

0.75* 

-0.72* 

-0.22 

^star 

0.67* 

-0.53* 

0.09 


Note. — Stars symbols indicate significant cor¬ 
relations (p-value < 2.5%) For a typical num¬ 
ber of 45 fully detected data points, p-values 
are < 50%, 10%, 5%, 2.5%, 1%, and 0.1% for 
\r^\ Z 0.10, 0.25, 0.30, 0.34, 0.39, and 0.49, re¬ 
spectively. 


trends across the ~ 2.1 dex (a factor of 125) spanned by 
both ISM parameters. As either Smoi/^iatom or Sstar increases, 
/dense increases by a factor of > 4 (fitted power-law indices 
~ 0.3; see equations in the top comers of the panels), whereas 
SFE^ense decreases by a factor of ~ 8 (power-law indices 
~ -0.4). In contrast, S appears virtually independent 

of Snioi/diatom Sstar within the scatter, in rough agre emen t 
with more detailed studies in HERACLES galaxies (|L13|l. 
These trends are confirmed in Table|^ where we list the corre¬ 
sponding rank coefficients. Clearly, Smoi/^iatom and Zstar show 
strong (anti)correlations with /dense and S FE^f^nse but no sig¬ 
nificant correlations with S FE^oi. 

Eig. S argues strongly that some intrinsic properties of 
molecul^ clouds vary systematically across galaxy disks. Eor 
our default assumption of fixed CO and HCN conversion fac¬ 
tors, we find that /dense and ^EEdense show opposite trends, 
in the sense that the S EEdense decreases and /dense increases 
as we move from low-Sstar. low Smoi/diatom disks towards high 
high Stnoi/Satom galaxy centers. Erom Equation 10 this 


renders S FE^oi almost uncorrelated with the local conditions 
considered here. In Appendix [A| we demonstrate that these 
results are robust against changes in the CO transition used to 
trace molecular gas (i.e., 7=1-0 instead of 7 = 2 - 1) or in 
the chosen SER tracer. 

Cast in observational terms, our results show that the ratio 
of HCN-to-CO intensity (oc Sdense/^moi) increases systemat¬ 
ically from a typical value of 1/30 (in K) in the disks of 
normal star-forming galaxies to a higher ratio of ^ 1 /10 in ac¬ 
tively star-forming galaxy centers and starbursts, i.e., a factor 
of ^ 3.5 increase over the two orders of magnitude in Sjtar and 
Smoi/^atom that we Study. At the same time, the IR-to-HCN 
ratio (oc SsFR/Sdense) drops, by a factor of =» 6 over the same 
range, so that active galaxy centers appear to display system¬ 
atically less infrared emission (thus, lower star formation rate) 
per unit dense gas emission. These two trends cancel, though 
not entirely, so that the ratio of IR-to-CO (oc ZspR/^imoi) varies 
more weakly over the same range, systematically changing by 
less than a factor of 2. 

Broadly, these results have two straightforward interpreta¬ 
tions: either S EEdense and /dense indeed change as described, 
or physical conditions within the molecular clouds change 
systematically in such a way that the conversions of CO and 
HCN intensity to total and dense molecular gas produce the 
trends that we observe. In either case, the observations im¬ 
ply that the physical conditions in the molecular ISM, and 
in molecular clouds specifically, change systematically across 
galaxy disks, but in a way that conspires to yield only weak 


apparent variations in the SER per unit molecular gas. We 
explore these two interpretations in the next sections. 

5. COMPARISON TO STAR FORMATION MODELS 

Because dense gas is believed to play a pivotal role in star 
formation, our results on /dense and ^EEdense have implica¬ 
tions for some of the most common models of star forma¬ 
tion applied to galaxies. In this section, we compare our ob¬ 
servations to exa mples of the two broad classes of current 
models: fi rst (Sect.|5.1|l, the s traightfo rwar d density-threshold 
model by Gao & Solomon| ( 20()4a|b i and Wu et al. <20051; 
second (Sect. |5.2|l, turbulent whol e cloud models ( [KMK05| 
KT07t Federrath & Klessen|2012| l. We specifically focus on 
the theor y o f turbule nce-regulated star formation discussed by 
KMK05|and|KT07| but refer the reader to a wider summary of 
whole-cloud mode ls by Federrath & Klessen ( 2012[ ). The al¬ 
ternative model by Narayanan et al.| ( |2008| l is based on hydro- 
dynamical simulations of entire galaxies, but shares a similar 


framework and results in similar predictions to the KMK05 

5.1. Density Threshold Models 

In a density threshold model, the star formation rate in 
a region depends chiefly on the available mass of molecu¬ 
lar gas above a certain volume density. The simplest form 
of such a model, and the most commonly adopted in extra- 
galactic research, assumes an environment-independent, fixed 
star formation efficiency of the dense gas traced by, e.g., the 
HCN(l-O) line (our ^EEdense)- Such a view olfers a sim¬ 
ple explanation for the almost linear correlation between the 
HCN(l-O ) and IR luminosi t ies of gal axies found in the early 
studies by|Gao & Solomon|(|2004a|bll. That relation contrasts 
with the (globally) superlinear relation between the CO(l-O) 
and IR luminosities when (U)LIRGs are included. Within this 
paradigm, the IR-molecular relations in galaxies emerge be¬ 
cause of a higher /dense in the IR-brighter objects. In sup- 


port of a density threshold model for star formation, Wu et al. 

(20051 concluded that galaxies and Galactic dense cores with 
Ltir > ^ Lq (as required to get reliable SER estimates 

from the IR continuum) align along the same HCN(l-0)-IR 
relation. Other studies of Galactic clouds based on dilferent 
obser vational methods have been found in qualitative agree¬ 
ment ( |Lada et ari|2012t [Heiderman et al.|[2010[ [Evans et al.| 
2014 ) 1 , favoring this simple picture in which the dense gas 
reservoir is the main determinant of the star formation rate. 

At face value, our results in galaxy disks are clearly at odds 
with these simple threshold models. We And that the S EEdense 
at kiloparsec scales systematically varies by a significant fac¬ 
tor of ss 6 across ~ 2 dex in S^oi/Satom Sstar- We do And 
that /dense Systematically increases by a similar factor moving 
from galaxy disks to galaxy centers, but the interplay of the 
two quantities means that S FE^o\ remains constant across the 
disks within the scatter. The apparently varying S FEisn^e in 
our observations means that the ability of dense gas to form 
stars must depend on environment and thus, presumably, on 
conditions in the host molecular cloud. This is at direct odds 
with the basic assumption of threshold models. 

One can envisage more complex density-threshold mod¬ 
els that could fit our resolved observations in galaxy disks. 
In principle, a model that allowed for variations in the den¬ 
sity threshold and/or in the star formation tim escale of the 
dense gas could solve the discrepancies (e.g., seejFederrath &| 
Klessen|[2012 l. However, making such adjustments quickly 
moves one into a gray area: if a density threshold depends 
sensitively on conditions in the parent molecular cloud and 
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Figure 2. Dependence on environment in our observations. From top to bottom, the dense gas fraction /dense (top row), the star formation efficiency of dense gas 
S -FEdense (middle row), and the star formation efficiency of molecular gas S FEmoi (bottom row) as a function of: the molecular-to-atomic ratio (Smoi/^atom. left- 
hand column) and the stellar mass surface density (Sstar^ right-hand column). The right-hand axis of each panel display the data in terms o f obs erved quantities. 
Filled and empty symbols correspond to detections and non-detections, respectively. The solid line shows our best fit to the data (Sect. \2A\ . Its equation is 
reported in the top right comer. 
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host galaxy, then it increasingly resembles the whole-cloud 
models. Alternatively, systematic variations in how CO and 
HCN trace the molecular and dense molecular medium could 
partly account for the observed trends, salvaging a more uni¬ 


versal density-threshold model. In Sect. 6.2 we explore 


whether our observations could be consistent with a simple 
density-threshold model given plausible variations in the con¬ 
version factors. 

5.2. The Turbulence Regulated Whole-Cloud Model 

“Whole cloud models” offer an alternative to the univer¬ 
sal density threshold view. Here we use this term to refer 
to models in which conditions throughout a molecular cloud 
ultimately affect the efficiency with which dense gas forms 
stars. These models have been contrasted with threshold mod¬ 
els in the literature several times, with recent Milky Way re¬ 
sults invok ed to support a simple density threshold ( |Evans | 
et al.|2014|l, while arguments using low resolution extragalac 


tic data used to argue for whole cloud models (e.g.,|Krumholz 
& Thompson|2013 i. Even considering only analytic models, a 


wide suite of whole clou d models exist in the literature ( |Eed-| 
[errath & Klessen||2012[ ). These almost all consider the star 
formation efficiency per free fall time, S FE^, as the fiducial 
quantity and so encode a fundam ental dependence on den¬ 
sity (e.g.,|McKee & Ostriker|2007||. SFEg further depends in 
most models on the virial parameter, Mach number, magnetic 
field, and turbulent character within the clouds. 


5.2.1. KMK05 Model Overview 

With such a wide set of free parameters, it is outside the 
scope of our paper to compare our data to the full set of 
models or possible cloud conditions. Indeed, because in our 
survey we have very limited knowledge of the cloud popula¬ 
tion at sub-beam scales, we lack many of the measurements 
needed for such a comparison. Instead, the main thrust of our 
treatment of whole cloud models is to ask: can such models 
plausibly create the kind of correlated variations in S EEdense 
and /dense that we see in our data? To do so, we focus on one 
of the earliest whole cloud models, th e theory o f turbulence- 
reg ulated s tar formation put forward by |KMK05] In a later pa¬ 
per, [KTOT] combined their theoretical treatment with a simple 
radiation transfer modej^that implicitly provides the conver¬ 
sion factors between gas masses and the observed molecular 
line luminosities (e.g., HCN(l-O)). To simplif y the comp ar- 
ison, we refer here to the original equations by KMK05 and 
treat con version fa ctors as fixed. 

In the |KMK05| theory, the star formation within a virial- 
ized molecular cloud with a lognormal distribution of den¬ 
sities depends primarily on the average gas density (h) and 
the intracloud turbulence. The latter is parameterized by the 
Mach number, M. The sense of the dependence is that h sets 
the characteristic timescale for star formation, the free-fall 
time, so that a high density implies more efficient star for- 
mation per absolute time interval. Meanwhile in the |KMK05 
model, a higher degree of intracloud turbulence, i.e., a higher 
M, tends to inhibit star formatiorp^ 

Eor our purposes, the KMK05| model boils down to the fol¬ 
lowing Equations (complemented by Equation [T^: 


Code available at http://www.ucolick.org/~krumholz/codes/ 

Note that this feature is not shared hy some other models, in which 
higher Mach numbers tend to drive higher rates of star formation rates; see 
[Federrath & Klessen|j2012|. 
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Here, rjf is the free-fall timescale evaluated at the average vol¬ 
ume density (h), which is proportional to 1 / Vh; esF is an effi¬ 
ciency parameter of the order of ~ 1.4% that depends on the 
virial state of the clouds; crpop is the width of the lognormal 
probability distribution function of gas density within a cloud, 
deter mined by th e Mach number {M). 

The |KMK05] model has two key features. Eirst, unlike 
a simple threshold model, whole cloud models predict both 
/dense and SFE^oh and by extension SFEis„se- Although 
these models have more free parameters, they also attempt 
to explain more of the star formation process; by contrast, 
the density threshold models push the topic of regulating star 
formation into the f ormation o f dense gas, which they do not 
address. Second, the |KMK05] model has two different mecha¬ 
nisms to alter /dense: the fraction of dense gas increases along 
with both the average density, h, and the turbulent Mach num¬ 
ber, h\. The latter dependence occurs because a higher Mach 
number broadens the density PDF (Eq.[T3|l. Consequently, for 
a fixed average density, it also increasesme mass of gas with 
densities above the cutoff for the HCN-traced gas, Wdense (this 
happens as long as ndense > «, wh ich is always ex pected for 
plausible conditions; see Sect.|2.2|l. In the |KMK05|model, the 
variations in h and Ad affect the star formation efficiencies in 
different ways, however: 

• If h increases, Tff decreases, enhancing S FE^aX- In con¬ 
trast, the product /dense X Tff changes little across a plau¬ 
sible range in h, which would render S EEdense almost 
constant. 

• If Ad increases, /dense increases, while SFE^aoX would 
moderately decrease. Combining these, ^EEdense 
would experience a more significant drop (due to the 
additional 1 //dense factor). 


Based on this model, |KT07| (see also Sect. |7|i posit that an 
increase by 2-3 orders of magnitude in h along the IR se¬ 
quence of galaxies could mostly explain the different power- 
law indices for the IR-CO(l-O) and IR-HCN(l-O) relations in 
galaxies. While a high Mach number in (U)LIR Gs com pared 
to normal SF galaxies could also come into play, |KT07 1 argue 
that this parameter would likely play a secondary role. 

5.2.2. Results 

Because we do not know the detailed structure of molecular 
clouds within our beams, the agreement between our data and 
Equations 1 11 1 13[ depends on some assumptions about four pa¬ 
rameters that are not strongly constrained a priori: the CO and 
HCN conversion factors, the density cutoff for HCN emission 
(nd ense)5 3^ri(l CSF- To ni3,tch thc rCgioil of S FEy^qX S ^£^dense 
/dense parameter space spanned by ou r obser vations, we set 
«dense to 2.8 X 10^ cm“^, as assumed by|KT07 but adopt a esF 
efficiency parameter that is 1/5 times the val ue prefe rred by 
them. The same mismatch in esF was noted by|GB12 
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Figure 3. Comparison between the S F£mol-/dense (left-hand panel) and S FFdense-/dense (right-hand panel) relations in our data and in the|KMK05|model. The 
top and right-hand axes of each panel display the data in terms of observed quantities. For the data, symbols are the same as in Fig.fT] The power-law fit equation 
(gray line) is shown in the top-right comer. The model predictions are represented by colored lines. To obtain the purple curves, labelled by the assumed M. in 
the left-hand panel, we fix the Mach number and vary the average gas density. To obtain the orange lines, labelled by the assumed h in the right-hand panel, we 
fix the average density and vary the Mach number. 


Figure displays our data in the S /’£’nioi-/dense (left) and 
S Fit dense 7dense (fight) planes, with power-law fits shown as 
grey lines. We overlay the curves predicted for realistic con¬ 
ditions when we vary either h (fixed M — 25, 50, and 100, 
from top to bottom; purple) or M (fixed h - 10, 10^, 10^, and 
10“* cm , from left to right; orange). Thus the purple and or- 


ange grid shows a plausible parameter space from the KMK05 
overlaid on our observations. 

Figure shows that our data are better described by the 
orange tracks obtained by varying the Mach number (Ai) 
while holding the average density constant at 100 cm“^. 
The mismatch between our data and the curves for constant 
Ai (purple) is particularly apparent in the right-hand panel, 
which shows that the spread in 5FFdense that we measure 
is hard to reproduce only by varying the average density of 
otherwise identical clouds. As a simple way to quantify this 
(dis)agreement, we fitted the data in each panel with curves 
predicted by the |KMK05| model, fixing either n or AI while 
leaving the other parameter free. In the S F£dense-/dense space, 
the best-fit curve for fixed Ai has a four times higher reduced- 
than the best-fit curve for fixed n (in the S FE^oi - /dense 
plane the fits have more comparable quality). 

Thus a simple comparison to the whole-cloud models of 


KMK05 shows that, with reasonable inputs (and after adjust- 
ing the esp parameter), the models span a range of S FFmoi- 
SFEdense-/dense parameter space that ca n match o ur data well. 
For the specific functional forms of the |KMK05 1 models, our 

of the mode 


data are well fit by an implementation of the models in which 
clouds are significantly more turbulent, but not much denser, 
near the galaxy centers than in the disks. If this were the case, 
the trends that we find across galaxy disks would be differ¬ 
ent from those between SF galaxies and (U)LIRG, which m- 
pear to be chiefly driven by an enhancement in h (see Sect.|^. 
However, we caution that the specific numerical dependen¬ 
cies of the |KMK05| model do differ from other, more re¬ 


cen t, analyti c models, including those that otherwise adopt 
the |KMK05| approach but treat the free fall time on multiple 
scales. Broadly, the wider parameter space spanned by whole 
cloud models can accommodate our data in a way that a sim¬ 
ple threshold cannot. 

6. THE EFFECT OF CONVERSION FACTORS 

So far, we have assumed fixed. Milky Way-like conversion 
factors for CO and HCN. Under this assumption, our data 
strongly favored whole-cloud models over a universal density 
threshold. An alternative hypothesis would be that the chang¬ 
ing physical conditions within our clouds mainly change the 
conversion between CO and total molecular gas mass or HCN 
and dense molecular gas mass. Though an exhaustive study 
of the conversion factor is well beyond the scope of this paper, 
here we test a suite of plausible conversion factor variations 
and investigate how they could alter the simple picture dis¬ 
cussed above. Specifically, we build a grid of plausible con¬ 
version factors for CO and HCN and, at every point of the 
grid, we check whether our data, combined with these con¬ 
version factors, exhibit basic agreement with either of our two 
models. 


6.1. Background 

Significant deviations in a^o from its Milky Way value 
have been predicted and observed in g alaxies, both at globa l 
and smaller scales (see the review by |Bolatto et al.||20131). 
Broadly, aco appears to exhibit two rnain behaviors! First, 
aco seems to increase with decreasing metallicity (or dust- 
to-gas ratio), as CO emission becomes confined to a smaller 
fraction of a molecular cloud. Second, aco tends to drop 
in extreme star forming environments (e.g., galaxy centers, 
(U)LIRGs, an d starbursts) where th e gas is more turbulent and 
more excited. |Bolatto et al.|(|20131l suggested that this second 
trend could be usefully parameterized as a decrease in aco 
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with increasing total mass surface density (Stot = ^stai + 2igas), 
although Slot may not be its ultimate physical driver. 

Compared with aco, the HCN conversion factor, cthcn, has 
been much more poorly characterized. There are, however, 
practical reasons to expect that the two conversion factors fol¬ 
low each other to some extent. First, ancN must be lower in 
some (U)LIRGs than in SF galaxies to ensure that /dense < 1 ■ 
More physically, if cthcn emerges from a turbulent, opaque 
medium — even if it makes up only a portion of a typical 
cloud — we expect similar dependences on the degree of tur¬ 
bulence and excitation. Second, the density above which the 
HCN(l-O) line is effectively emitted could vary as certain 
cloud properties change (e.g., average density, temperature, 
turbulent velocity gradient). This would also make cthcn, that 
converts the line lumi nosit y into gas mass above a hxed den¬ 
sity cutoff, vary (Sect. [23. 

We adopt a simple, ad hoc approach to gauge the interplay 
of conversion factors, models of star formation, and our ob¬ 
servations. To do so, we simply assume that, to some very 
coarse approximation, both aco and othcn may vary as some 
power law function of Sstar, that is: 


ttHCNC^^star, THCn) - 


\ 100 Mo pc- 

^star 




\yco 


(14) 

(15) 


Here, the prefactors are set to the default Galaxy-like HCN 


and CO conversion factors defined in Sect. 2.2 Equations 
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and 15 are remini scent of the empirical law proposed by Bo- 
|latto et aT] ( |2013| l to match a set of aco v alues culled from 
the lite rature, including those estimated by Sandstrom et al. 
( 12013 1 across the disks of HERACLES galaxies. [Bolatto et al. 
( |2013 1 specifically proposed yco = -0.5 above a total mass 
surface density of 100 Mq pc'^ in addition to an exponen¬ 
tial metallicity dependent term. We neglect the metallicity 
term because our survey focuses heavily on the metal-rich 
parts of the HERACLES sample (we checked that this term 
would hardly affect our results if included). We approximate 
Stot ~ 2istar because, in our sample, the total surface density 
is almost always dominated by Sstar (usually by a factor of 
> 3). Einal ly, we have formally dropped the Stot threshold ad¬ 
vocated by|Bolatto et al.|(|2013|l because essentially all of our 
detections have Xstar above this value (see Eigure]^. 

We emphasize that Equations and [^represent a coarse 
approach to the variability of con version factors. The det ailed 
study of HERACLES galaxies by|Sandstrom et al.|(|2013|l sug¬ 
gests that aco m an individual galaxy can be better described 
by a flat profile over most of the disk and an abrupt drop near 
its center and we do not expect disk surface density to be the 
only (or even main) physical driver of conversion factor vari¬ 
ations. Nevertheless, power laws offer a useful tool for the 
kind of exploratory calculations that we carry out in this sec¬ 
tion. These are already highly approximate due to the lack 
of observational and theoretical constraints on cthcn and a 
power law formalism allows us to treat ancN and aco m a 
similar way with only a minimal set of free parameters (yco 
and 7hcn)- Moreover, we are specifically interested in how 
conversion factor variations could “rescue” a density thresh¬ 
old model. Eig.j^shows that, at least in the specific case of our 
set of pointings, ancN should roughly behave as a power-law 
function of Sstai in order to cancel the systematic gradients in 
IR-to-HCN that we observe. 


We can consider a grid of plausible yco yncN, with 
“plausible” defined according to the following rules: 

i) Most evidence in the literature point to lower conver¬ 
sion factors in regions of higher surface density (galaxy 
centers and (U)LIRGs), so we assume yco ^ 0. 

ii) We expect the fraction of gas in the molecular phase, 
traced by Smoi/diatom, to increase with the depth of the 
gravitational potential well, traced by Sstar- We And that 
we require yco > -0.82 to ensure that this ratio and 
Sstar are positively and significantly correlated and we 
consider only values above this value. 

iii) We expect that /dense increases along with the HCN/CO 
line ratio. We And this happens when yncN - yco > 
-0.42. That is, the Q’hcn/Q'co ratio cannot be so low 
at galaxy centers that we actually have lower dense gas 
fractions for higher HCN/CO ratios. 

iv) We require that less than 5% of our positions have 

a nominal /dense > / which is unphysical (even 

widespread /dense ~ 1 is implausible). Eor this, we And 
that the condition yncN - Tco < 0.58 is necessary. That 
is, the ancN/ffco ratio cannot be so high at galaxy cen¬ 
ters that all of the molecular gas is in the dense phase. 

The four conditions stated above deflne a rectangular region 
in the (yco, yncN-yco) plane that encapsulates plausible vari¬ 
ations in the CO and HCN conversion factors. The changes in 
yco and yncN will modify the mutual relations between /dense, 
SFE^oi, and inferred from our observations. Thus, 

exploring this parameter space, we study the compatibility be¬ 
tween our observations and models as a function of possible 



relevant for the comparison between observed and predicted 
trends that we carry out below. 


6.2. Density Threshold Models 

The key assumption of the simple densi t y-thres hold hy- 
pothesis supported by, e.g.,|Gao & Solomon](2004b|l andfWu 
|et al.| ( |2005| l is that SFE^ense is Axed within the scatter. In 
Figure !^ the orange-colored area within the thick black con¬ 
tour is me region of the (yco,THCN - yco) parameter space 
where this could be fulfllled: inside it, S FEdense (derived from 
our data and the conversion factors at each point) would not 
signiflcantly correlate with Sstar, nor with S^oi/^atom- Essen¬ 
tially, this region stretches along the straight lines for which 
-0.55 < yncN ^ -0.35, which is the cthcn required to can¬ 
cel out the gradients in S FEdense found in Sect, We remark 
that, to cancel these gradients, ancN should rou^ly behave as 
such a power law when restricted to our sample of pointings, 
regardless its true functional form. 

In more detail, the region without systematic variations in 
S FEdense consists of three different sections where the thresh¬ 
old model has a decreasing predictive power for our sample 
of pointings. In the flrst section, labeled “A” in the flgure, 
S FEmoi and /dense are mutually correlated and show system¬ 
atic variations with respect to either Smoi/^atom or ^star- Thus, 
a priori, systematic increases in S FEmoi could be ascribed to 
increases in /dense, as expected for a threshold model. Section 
“A” corresponds to values of -0.7 < yco ^ -0.3, i.e., com - 
patible with the yco = -0.5 value in Bolatto et al. (20131, 
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Figure 4. Compatibility between our data and the density-threshold models 
across the (ycOT yHCN ” Jco) region defined by our boundary conditions (see 
text; we only represent yncN ” Tco < 0-35 to keep the plotting box square). 
Within the orange area, the observations and the density-threshold model ai‘e 
compatible {S -FFdense is neither significantly con'elated with Sstai- nor with 
^moi/^atom^ see text for details). The purple contours enclose the regime there 
are systematic variations in either and /dense. but both variables are 

either anticorrelated or not significantly correlated at all. The purple crosses 
indicate the intersection between this regime and the region where S FFdense 
shows no systematic variations. This defines section ”C”, while sections ”A” 
and ”B” are located at either side of it. 


and to variations in ancN that are largely in lockstep with 
those in aco (Ithcn - Tcol ^ 0.2), perhaps in response to 
the same physical mechanisms. The significant correlation 
between /dense and S FEi„o\ arises because the Sstar dependent 
term in othcn flattens ^FiSdense. whereas the steep trend in 
/dense is barely affected as variations in (Thcn and aco cancel 
each other out. 

In a second section, labeled “B”, the systematic gradients 
in /dense. S Fiidense. and SFE^oi across galaxy disks all cancel 
out. This leaves the density threshold theory plausible, but 
essentially unprobed within the scatter of our data. Section 
“B” corresponds to only weak variations in aco across galaxy 
disks, but much stronger gradients in ancN- Essentially, anew 
varies in such a way that it cancels the observed variations in 
the IR-to-HCN and HCN-to-CO ratios, whereas the trend in 
the IR-to-CO ratio remains unchanged. 

Finally, there is a third region, labelled “C”, where there 
are systematic variations in either SFE^^i and /dense, but both 
variables are either anticorrelated or not significantly corre¬ 
lated at all. In this case, /dense could not be considered a major 
driver of SFE^oi, which would rather be regulated by an hy- 
potethical hidden variable. Therefore, even though S FE^ense 
would remain constant, the density threshold model would not 
be not predictive of SFE^^i. 


6.3. The Turbulence-Regulated Model 

To illu strate the impact of changing 7 hcn and yco on tbe 
|KMK05| model, we consider two extreme cases: one in which 
we keep h fixed while v arying M and one in which we fix 
M. and vary h (see Sect. 5.2.1 and Fig. |^for an illustration 


of such models). Of course, in a real systemboth quanti¬ 
ties may vary and, as discussed above, the |KMK05| model 
represents only one implementation of a turbulence-regulated 
whole-cloud model. The essential features of the models, 
which our data should reproduce to be compatible with the 
theory in a region of (yco. Thcn - Tco) space, are: 


• For the fixed M case, both /dense and S FE^o\ should in¬ 
crease when h increases; that is, we require a significant 
positive correlation between /dense and S FEmo\- 

• For the fixed h case, as Figure shows, /dense increases 
when A4 increases, whereas SFE^oi decreases more 
mildly; as a result, we expect either an anticorrelation 
between /dense and S FE or even no significant cor¬ 
relation given the observational uncertainties. 


We implement these qualitative behaviors as requirements 
on the p-value of the rank correlation coefficient between 
/dense and SFE^oi and show the valid regions in (yco. Thcn - 
yco) space in Fig.|^ Fig.j^shows that almost all of the plau¬ 
sible (yco. Thcn) parameter space can be covered by some 
variant of the turbulence-regulated model. Under the assump¬ 
tion of fixed h (left-hand panel), our observations agree with 
this theory when yco ^ -0.2 or yncN - Too ^ -0.2. This in¬ 
cludes the case of fixed conversion factors (yco = Thcn = 0 ) 
adopted throughout the first part of the paper. If M remains 
fixed while h varies (right-hand panel), observations and mod¬ 
els agree within the complementary region, but also in a small 
patch at yco - 0 and yncN - Too - -0.4. This is the patch 
where the CO and HCN conversion factors conspire to elim¬ 
inate all gradients (/dense. SFE^oi, and ^FFdense) from our 
sample. Of course, inasmuch as this conspiracy yields a ho¬ 
mogenous population of star-forming clouds throughout the 
sample, our observations could be compatible with any cloud- 
scale star formation theory. 


6.4. Conclusions 

Our exploration of plausible aco and ancN shows that com¬ 
patibility with a simple density-threshold model of star forma¬ 
tion places specific, non-trivial requirements on the behavior 
of the CO and HCN conversion factors. Essentially, ancN 
must vary in a way that cancels out the observed variation 
in the IR-to-HCN ratio and its value must be fairly tightly 
coupled to aco- Fy contr ast, the turbulence-regulated whole 
cloud model of |KMK05|can be compatible with almost any 
region of aco-o'HCN space that we explore. To some degree, 
this reflects the larger suite of free parameters in the whole 
cloud models, but this does not entirely invalidate the con¬ 
trast: it may very well be that more physics, and so more free 
parameters, are needed to explain the observations. 


7. DISCUSSION 

We interpret our results to offer a challenge to simple den¬ 
sity threshold models, thoug h not to c onclusively rule them 
out. Whole-cloud models like IKMK05 l have more success ex¬ 
plaining our observations because they can span a wider part 
of 5T'£’[noi-/dense-5F£’dense space. We discuss shortly the im¬ 
plications of combining each of thes e mode ls with our data 
and the unresolved observations from lGBO 


7.1. Implications for Density Threshold Models 
Density-threshold models found support in the remark- 


able linear IR-HCN(l-O) relation found by Gao & Solomon 
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fixed n, M. varies 


fixed A1, n varies 



- 0.8 - 0.6 - 0.4 - 0.2 0.0 

Yco 



Figure 5. Compatibility between our data and the |KMK05| model for variation in the Mach number (left-hand panel) and in the average density (right-hand 
panel). Each panel covers the (yco. THCN “ yco) region dehned by our boundary conditions. The orange area indicates the (yco. yilCN “ yco) values for which 
observations and theory are compatible (see text for details). 

threshold; the simple correlation of HCN and IR luminosity. 
Moreover, an ad hoc tuning of anew cannot address the higher 
apparent S FfSdense in (U)LIRGs, since the basic expectation is 
that their ancN should be lower than that of SF galaxies. Ap¬ 
plying such a correction would exacerbate, rather than allevi¬ 
ate, the differe nce in S FEjense between both kinds of galax ies 
(for details, see GB12| l; but see Papadopoulos et al. ( 2012|l for 
an alternative discussion of conversion factors in (U jLlRGs. 


( 2004a| l and Wu et al. ( |2005[ ), which encompassed systems 
covering ~ 7-8 dex in Ltir- This correlation represents al¬ 
most the only extragalactic observational support for such a 
model. In plots that span such an enormous range in lumi¬ 
nosity, our data do broadly “fill in” the previously missing 
part of this correlation. However, in detail our data also ex¬ 
hibit clear, systematic variations in the IR-to-HCN luminos¬ 
ity ratio. That is, luminosity-luminosity plots or power law 
fits spanning many decades would miss significant, system¬ 
atic environmental trends that appear as “scatter” around the 
luminosity-luminosity scal ing. We suspect that along with in¬ 
trinsic cloud-cloud scatter ( |Ma et al.|2013| l and observational 
uncertainties, the wide range in measured SFR/Lucn ratios 
(> 2 dex and ~ 1 .4 dex in Galact ic cores and galaxy aver¬ 
ages, respectively; Wu et al.|2005| GB12) may be driven by 
the same systematic variations that we see in this survey. 

Other studies have highlighted variable SFE^ense' recent 
studies at whole-galaxy scales have pointed out that the aver¬ 
age SFEiense of (U) LIRGs is at least 3-4 times higher than 
that of SF galaxies (|GB12| see also |Gao et al.||2007|l, un¬ 
like what |Gao & Solomon| ( |2004a| l concluded. Our own fit 
to the GB12 data implies that the average 5 FT;dense would 
increase as Moreover, studying t he center of our 

own Milky Way galaxy using ammonia lines, [Longmore et al.| 


( |2013 1 demonstrated this high surface density region of a nor¬ 
mal disk galaxy to show a systematically low S FE^^„se, the 
same sense of the trend (^TEdense anticorrelates with Sstar) 
that we observe in this paper. 

For a fixed HCN-to-dense gas conversion factor, these re¬ 
sults are clearly at odds with the assumption of a fixed 
SFE^ense that Underlies threshold models. We show that 
our resolved observations could be compatible with a fixed 
(within the scatter) S Ffidense if aco and <thcn varied in some 
plausible way. However, we emphasize that doing so would 
remove the major extragalactic argument for a simple density 


7.2. Turbulence-Regulated Models and a Simple 
Observational View 

Whole cloud, turbulence-regulated models can produce a 
wide range of points in the 5T £’moi-./dense-^ F'fi’dense space. 
The specific model of |KMK05 1 therefore appears able to re¬ 
produce our data for almost all assumpt ions about conversion 
factors, though recall that, as in |GB12| we had to rescale the 
absolute scale of the star formation rate in the model (i.e., 
the efficiency parameter esp)- We saw that, for the case of 
fixed conversion factors, our data ag ree particularly well with 
an implementation of the KMKOS] model in which the aver¬ 


age cloud density remains approximately fixed while the tur¬ 
bulent Mach number increases with increasing stellar surface 
density. 

This primary dependence on the Mach number may seem 
at odds with the normal explanation for increased S FE^o\ in 
galaxy centers, which would attribute the contrast between 
galaxy disks and (UjLlRGs mainly to gas density. We add ress 
this in Fig. where we combine our data with the GB12| data 
in the mostly observational spaces 2sfr//co vs. /ncNio/fco 
(left-hand panel) and 2sfr//hcnio vs. /ncNio/fco (right-hand 
panel), where Ssfr is estimated linearly from the IR. This fig¬ 
ure encapsulates our best current observational understanding 
of HCN, CO, and star formation in galaxies. In both panels, 
the combined data sets span a contiguous region of the param- 
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Ihcn/Ico 



Ihcn/Ico 


Figure 6. Comparison between observations and the|KMK05[model : Ssfr/^co (left-hand panel) and Xsfr/^hcnio (right-hand panel) as functions of /hcnio/^co 
in our data and the GB12 extragalactic sample (note that intensity and luminosity ratios are formally identical, which allows us to include in this figure those 
galaxies of the GB12 sample that lack an estimate of the molecular disk area). Symbols as in Fig.^ For the sake of visualization, limits are not displayed. The 
purple lines correspond to lines of constant Mach number At = 25, 50 and 100, as indicated by the Tabels of the same color. They are obtained by increasing the 
average gas density (n) from 10 cm“^ to 10^ cm“^ in steps of 0.5 dex. The gray lines correspond to constant n conditions, as indicated by the brown labels (for 
clarity, only one in two lines is labeled). The space between them is colored in green-orange-red hues, with reddish ones corresponding to higher densities. To 
convert (dense) gas masses into line luminosities, we assume the default HCN and CO conversion factors of this paper. The gray area in the background indicates 
the range of observables that can be reproduced by the|KT07|model, where non-fixed HCN and CO conversion factors are derived from a simple radiative transfer 
code coupled to the |KMK05| cloud model. The upper and lower envelopes of this area correspond to the normal and starburst set of parameters defined by |KTQ7| 
They assume Mach numbers of 30 and 80, respectively. 


eter spac^^ Remarkably, there is no obvious separation be¬ 
tween (U)LIRGs and SF galaxies (i.e., no hint of bimodality), 
unlike what Fig. [T| suggeste d. Fig. [6|shows that, for SF galax¬ 
ies in both our sample and |GB12| there is no strong correla¬ 
tion between Zsfr//hcnio and /hcnio/7co (right-hand panel). 
Such a correlation exists, but only in the contrast between SF 
galaxies and (U)LIRGs and more weakly within the (U)LIRG 
sample . At the same time, both sets of SF galaxies (ours and 
GB12|l show a wide range in Ssfr//hcnio. including an anti- 
correlation with /hcnio/^co that becomes much clearer and 
unambigous when we consider the Sstar as the independent 
variable in our own sample. 

Colored region s in the ba ckground compare all data to the 
predictions of the |KMK05] model. These are the same shown 
in Fig. except that predicted surface densities of (dense) 
molecu& are converted to line intensities via our default hxed 
conversion factors. The regions with constant h are colored 
from green to red, with red hues indicating higher densities, 
for the sake of visualization. 

In the /hcnio/7co vs. Esfr/7co plot (left-hand panel), the 
collection of data from SF disks through (U)LIRGs tend to 
follow the purple lines (constant Mach number), with most 
but not all of them spreading between AI = 25 and Ai - 
100. Therefore, as we move from SF galaxies to (U)LIRGs, 
7hcnio/7co and 2sfr/7co would mostly increase because h 
changes from 100 cm“^ (solid line near the bottom-left cor- 

The only outlier is APM 08279+5255, a luminous BAL quasar at red- 
shift z ^ 3.9 with an exceptionally high Lhcn/^co ratio. Multitransition 
studies of this source has shown that its HCN excitation is driven by IR pump¬ 
ing to a significant extent ^Weil3 et al.|2007), in which case the line luminosity 
is not a reliable tracer of gas mass. 


ner) to ~ 10"^ cm“^. In this plot, changes in Ai would be sec¬ 
ondary and would mostly introduce scatter around the main 
trend. 

The /hcnio/7co vs. Esfr/7hcnio plot in right-hand panel 
provides us with a complementary view. Here, the most rele¬ 
vant trend comes from considering only our resolved obser¬ 
vations of SF disks (gray dots), which spread out roughly 
along lines of hxed n, mostly between « = 10' ^ cm“^ and 
h = 10^^ cm“^ (green-yellowish areas). Therefore, as we 
move across galaxy disks towards their centers, /hcnio/7co 
would increase and 2sfr/7hcnio would decrease mostly be¬ 
cause Ai increases from 25 to > 50. Changes in n would 
play a secondary role here. A less prominent but signihcant 
trend is found when SF galaxies and (U)LIRGs are globally 
compared. In this case, the average difference in S FE^ense be¬ 
tween the two families would mostly stem from a difference 
in the Mac h number. 

Thus the |KMK05 model provides a qualitative picture that 
can reasonably account for the observations. We caution that 
the quantitative details of this picture are far from settled, 
however. Firstly, we had to rescale their model to produce 
the comparison grids, and the exact numbers remain imper¬ 
fect matches to typical average cloud properties; e.g., typi¬ 
cal Milky Way clouds have AI ^ 10 and a volu me averaged 
density h ~ 100 cm~^ (e.g., Heyer et al.||2009 i, not perfect 
matches to the models. Secondly, a signihcant fraction of data 
points lie outside the generous 25 - 100 range in Mach num¬ 
bers that we have assumed to generate the comparison grids. 
Finally, the model requires the clouds to be less turbulent in 
(U)LIRGs than in SF galaxies, which is not the sense seen 
in observations. Still the sense of the trends seems reason- 
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able; several studies suggest that clouds in galaxy centers are 
dens er but also more turbulent than those in the g alactic disks 
(e.g., Oka et al.|2001[ Rosolowsky & Blitz|2005 Leroy et al., 

submitted). _ 

These inconsistencies between the IKMK05I model and the 
extragalactic data might stem from our crude assumption of 
fixed CO and HCN conversion factors for all data points. 
They might be partly alleviated by solving the radiative trans¬ 
fer independently in each pointing/source to calculate more 
accurate conversion factors. To illustrate this, we plot in 
Fig. |6l as a grey background the predictions from the |KT07 
modd, which combines the |KMK05| model with a basic ra¬ 
diative transfer model. The gr ay area is limited by two sets 
of cloud conditions defined by |KT07| one tries to reproduce 
typical conditions in SF galaxies, and assumes, in particular, 
Al = 30 (upper envelope); the second represents conditions in 
starburst galaxies, and assumes Al = 80 (lower envelope). In 
both cases the average gas density (h) is left as the only free 
parameter. We globally shifted the predictions in the vertical 
direction to better match the data clouds. The most apparent 
result is that, once the radiative transfer is taken into account, 
the model can explain the observed range in £sfr/7co and 
2isFR/7HCNio ratios with a narrower range of Mach numbers 
(30-80). 


In spite of this relative improvement, the KT07 model nei¬ 
ther reproduces the data to a satisfactory extent, though. E.g., 
it still requires a too broad range of Mach numbers and it pre¬ 
dicts clouds in (U)LIRGs to be less turbulent that those in SF 
galaxies. These two flaws might stem from th e simplifica¬ 
tions in the radiative transfer code use by KT07 (see the paper 


for details). They may also be solved by combining radia¬ 
tive transfer models with the more recent turbulent models in 


Federrath & Klessen (20121 and following. In general, the 
assumption of virialized clouds with lognormal density PDF 


implicit in the KMK05 model does not necessarily hold in the 
disturbed ISM of many (U)LIRGs. 

7.3. Simple Predictions and Logical Next Steps 

The contrast of our data with both models indicates logical 
next steps. First, while the CO-to-H 2 conversion factor has re¬ 
ceived substantial attention, the conversion between HCN and 
the dense gas mass remains far more uncertain. Especially the 
dependence of this value on environment is very uncertain. 
We expect HCN (1-0) and similar high critical density, low 
excitation lines (e.g., HCO+ (1-0), CS (2-1), or HNC(l-O)) to 
be the most accessible tracers of dense gas in external galax¬ 
ies for the foreseeable future. A quantitative understanding of 
how qthcn and similar quantities depend on environment will 
be essential. New observations, e.g., of optically thin isotopo- 
logu es over a wide range of env ironments, are clearly needed 
(e.g., Papadopoulos et al.|2012 1 . 

Along similar lines, the success of the “whole cloud” mod¬ 
els hinges on the properties of individual turbulent molecu¬ 
lar clouds. These properties, including average density and 
Mach number, can be measured and compared among envi¬ 
ronments. In the same way that knowing cthcn will provide a 
rigorous constraint on threshold models, folding cloud prop¬ 
erties into the comparison to whole-cloud models will vastly 
constrain the area of viable parameter space. A number of cur¬ 
rent mm-wave facilities are able to resolve i ndividual clouds 

at the distances of many of our targets (e.g. IColombo et al.l 

- 

8. SUMMARY 


We have used the IRAM 30-m telescope to observe the 
HCN(l-O) line in 62 positions (48 detections) across the 
disks of 29 nearby galaxies also targeted by the HERA¬ 
CLES CO(2-l) survey. Combining these data with previ¬ 
ous multiwavelength observations, we study the relations be¬ 
tween the dense molecular gas traced by the HCN emission 
(n > 10^ - 10^ cm“^), the bulk of the molecular gas (traced 
by the CO data), and the star formation rate (SER; traced by 
infrared emission). We highlight that; 

• Unlike most previous studies of HCN in galaxies, our 
observations span a wide range of galactocentric radii 
(from galaxy centers up to ~ 80% of the optical radius) 
at linear resolution (~ 1.5 kpc on average) that isolates 
local properties of a galaxy disk while still averaging 
over a large population of clouds. This allows us to 
probe a wide range in environmental parameters. 

• We compare our data set to a compilation of unre¬ 
solved observatio ns of s tar-forming (SE) galaxies and 
(U)LIRGs (from |GB12| l. Our resolved observations 
seamlessly expand the range of surface densities sam¬ 
pled by HCN studies of extreme systems to lower val¬ 
ues. They also fill in an intermediate, so far largely un¬ 
explored, luminosit y regime in the HCN-IR luminosity 
correlation found by |Gao & Solomon| ( |2004b| l; |Wu et ah] 

We focus our analysis on three parameters; the dense gas 
fraction (/dense HCN/CO), the star formation rate per unit of 
dense molecular gas (S FEdense SER/HCN) and the star for¬ 
mation rate per unit of molecular gas (5F£dense SER/CO). 
The last two are often referred to as the “star formation effi¬ 
ciency” of the dense and total molecular gas. We study how 
these quantities relate to conditions in the disk. We focus on 
the molecular-to-atomic mass ratio (Emoi/Eatom) and the stel¬ 
lar mass surface density (Estar) as useful indicators of the lo¬ 
cal conditions in the disk. Both quantities are highly covariate 
with galactocentric radius. We find that; 

• Eor the simplest assumption of fixed luminosity-to- 
mass conversion factors for CO and HCN, we find 
that /dense HCN/CO) tightly correlates with both 
Smoi/Satom and Egtar (Spearman’s rank coefficients 0.67 
and 0.75). Across ~ 2.1 dex in both disk structure pa¬ 
rameters, /dense Systematically varies by a factor of ~ 4. 

• In contrast, 5F£den.se IR/HCN anf/-correlates 
with Emoi/Eatom and Sstar (rank coefficients -0.72 and 
-0.53), systematically decreasing by a factor of ~ 6 - 8 
as we move from disk positions towards galaxy centers. 
The two trends cancel each other out, so that S FE^o\ 
shows no systematic variation with either Emoi/Satom 
nor with Estai. This implies that, to some degree, the 
apparent constancy of SFE^o\ known from previous 
studies arises from a conspiracy, as the sub-resolution 
cloud populations must vary systematically to produce 
the changing HCN-to-IR and HCN-to-CO ratios. 

We compare our results to two families of standard models 
of star formation in galaxies; 

• Eirst, we consider a simple density-threshold model, 
which, in its simplest version, assumes that the SER 
is proportional only to the mass of dense gas above a 
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certain density threshold traced by the HCN emission 
(thus, a fixed ^FEdense)- Such studies have been mo¬ 
tivated by local studies of Galactic clouds and the ap¬ 
parent linear correlation of HCN and IR across many 
decades in luminosity. Such a model is clearly at im¬ 
mediate odds with the apparent systematic variation in 
SFEdense across galaxy disks. Further, it predicts a 
strong, direct dependence of S FE^oi on /dense that our 
data do not immediately show. For fixed conversion 
factors, our data do not appear compatible with a sim¬ 
ple threshold model and contradict one of the main ob¬ 
servational pieces of evidence used to argue for such 
data. 

• In order for our observations to agree with a simple den¬ 
sity threshold model, the HCN-to-dense gas mass con¬ 
version factor, cthcn, must vary in such a way as to can¬ 
cel out the observed variations in the HCN-to-IR ratio. 
We explore the topic more quantitatively and show that 
only a very restricted range of possible combinations 
of CO and HCN conversion factors can satisfy such a 
model. 


We contrast the density threshold model with “whole 
cloud models,” focusing on the theory of turbul ence- 
regulated star formation formulated by KMK05I as an 


example (see also KT07 


and, for a wider view of such 


models, Federrath & Klessen||2012| l. In the |KMK05| 
model, SFE^oi increases with the average density of 
the molecular gas (n), whereas it is inhibited to some 
degree by increased turbulence (high Mach number, 
M). Meanwhile, /dense increases with both n andy^ 
We find that, for fixed conversion factors, the |KMK05| 
model can reproduce our data well, provided that the 
increase in /dense towards galaxy centers is driven by 
variations in Ai rather than in n. This good agreement 
exists for almost the full range of plausible conversion 
factors that we consider, reflecting the flexibility of the 
models. 


We interpret these results to represent a challenge to sim¬ 
ple density threshold models, though they do not rule them 
out. The simple correlation between IR and HCN luminos¬ 
ity spanning from starbursts to Galactic cores has represented 
an important piece of support for such models, our survey of 
galaxy disks shows that the main sites of star formation in the 
local universe do not obey such a simple relation (see also a 
dedicated study of M 51 by Bigiel et al, in prep.). 

The appendices report more details of the observations and 
investigate the effects of varying how we trace the SFR and 
molecular gas on our conclusions, showing that they remain 
robust to substitution of CO transition or SFR tracer. 
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Table A1 

Rank coefficients for /dense. S FEicme and S -F£mol as a function of local 
parameters, when either the CO(2-l) or the CO(l-O) transitions are used to 
estimate Zmoi 



Emoi from CO(2-I) 

5^mol from CO(l-O) 

/dense 

S REdense 

S FE^aX 

/dense 

S REdense 

S REniol 


- 0 . 59 * 

0.43* 

-0.07 

- 0 . 55 * 

0.43* 

-0.10 

^mol/^atom 

0 . 75 * 

-0.72* 

-0.22 

0 . 77 * 

-0.68* 

0.00 

^star 

0.67* 

-0.53* 

0.09 

0.76* 

-0.53* 

0.34* 


Note. — Star symbols indicate significant correlations. 
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APPENDIX 


A. ALTERNATIVE CO AND SFR TRACERS 
A.l. Checking 

By default in this paper, we assume a fixed CO(2-l)-to-CO(l-0) line ratio f? 2 i = 0.7 to estimate S^oi directly from the 
intensity of the CO( 2-l) line (Equatio n!^. There is evidence that /? 2 i could systematically vary across galaxy disks to some 
extent, however (see |Leroy et al.|2009| for the case of HERACLES galaxies; also references therein). In particular, R 21 could be 
up to ~ 1.5 times higher near sonie galaxy centers, as the gas excitation tends to be more efficient. By neglecting these var iations , 
we could be overestimating Smoi near galaxy centers. Thanks to our complementary observations of CO(l-O) (Sect. |2.1.1[ ), 
available for most of the observed positions, we can confirm that the assumption of a hxed R 21 has a minimal effect on our 
results. To do this, we derive an alternative 2n,oi after substituting R 21 = 0.7 in Equation 1^ with the line ratio derived from the 
CO(l-O) and CO(2-l) observations in each observed position. The CO(2-l) spectral cubes were previously convolved to the 
angular resolution of the CO(l-O) observations (22") to account for aperture effects. This implicitly assumes that the R 21 ratio is 
simila r at 22 " and 28" resolution. 

Fig. |Al| shows /dense. ‘SEEdense and SFE^oX as a function of Smoi/Satom when Smoi is derived from CO(l-O) (grey dots) and 
CO(2-l) (blue crosses). The plots where the variables are plotted against Sjtar are similar to these, so we do not show them here. 
In each panel, the cloud of blue crosses was slightly adjusted so that the two data clouds have the same med ian along both axes. 
This allows us to better compare the trends that we find for both CO transitions. It is apparent in Fig. |A 11 that those trends are 
nearly identical in each case. The fitted power laws are the same as those in Fig. [^within the errors, although they tend to be 
slightly steeper when Smoi is deri ved f rom the CO(l-O) intensity (i.e., by < 0.1). Tnis is consistent with a small increase in R 21 
towards the galaxy centers. Table |Al] shows that the rank coefficients discussed in Sect.fflare marginally sensitive to the chosen 
CO transition. The only noteworthy difference is that the correlation between SFE^^i ana Sstar is more significant for CO(l-O). 


A.2. Checking SspR 



Ssfr(L12) 
Mq Myr^* pc“^ 


= 634 




erg s“* sr ^ cm ^ 


Tcirrus-corrected 
^ 94 

COS(0 + 1.3-r COS(0. 

400 MJy sr-i 


(Al) 


In Equation |A 1 1 the measured 24 /rm intensity is corrected for the contribution from dust heated by pervasive radiation fields 
(’’cirrus”) unrelated to recent star formation. The cirrus component is proportional to the mass surface density of dust and a 
certain emissivity. Both parameters are derived from the dust models men tion ed in Sect. |2.2| A detailed comparison of different 
SFR tracers in the context of the HERACLES survey can be found in|L12|and|L13| 

We plot in Fig. |A2|the two star formation efficiencies against the molecular-to-atomic ratio, but now using the three alternative 
SFR tracers: the 24 jim continuum (left-hand column), the hybrid |L 1 2| calibration (middle column), and the Ha line (right-hand 
column). The plots against Sstar. not shown here, are similar to these. For the sake of comparison, we represent the results for 
our default SFR tracer (TIR) by means of blue crosses. These data points have been vertically shifted in each panel so that their 



































21 




Figure Al. /dense (left-hand panel), S FFdense (middle panel) and S FEy^a\ (right-hand panel) as a function of the molecular-to-atomic ratio. For grey dots, Xmoi 
is derived from CO(l-O) data as indicated in the text. For blue crosses, it is derived from the FiERACLES CO(2-l) spectra. 
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Figure A2. S FFdense (top row) and S FE^^i (bottom row) as a fu nction of the molecular-to-atomic ratio. From left to right, the SFR is derived from the 24 jim 
emission, a linear combination of 24 fim and Fin intensities jL12), and the FFcr emission. 


median Esfr is the same as that of the grey dots. This allows us to focus on trends, letting aside the relative calibration between 
the different tracers. 

The SFR calibrators have different levels of intrinsic scatter and are sensitive to different components of star formation. Some 
mild changes in the power-law fi ts can be expected, since the fraction of unobscured (i.e., Hcr-traced) star formation is known to 
increase towards fainter regions ( L12| l. In spite of all this, blue crosses and grey dots show similar trends (e.g., as captured by 
the fitted power-law indices) in most cases. We find only two noteworthy features, in the 5/’£’nioi(L12) and in the 5/’£’moi(HQ') 
panels: 


A few data points increase the dispersion in 5/’£’njoi(L12) at relatively low Smoi/Satom- On the one hand, there are three 
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Table A2 

Rank coefficients for S F^dense and S FE^oi as a function of local parameters, when different SFR tracers 

are used to estimate XsFR 



TIR (default) 

24/im 

L12 

Hff 

S F£”{16086 

^ ££mol ^ F £dense 

SfEinol ST^fidense 

‘^££mol ‘^££dense 

S F£:,^oi 

r25 

0.43* 

-0.07 

0.36* 

-0.05 

0.35* 

-0.13 

0 .68* 

0.34* 

^mol/^atom 

-0.72* 

-0.22 

-0.56* 

-0.08 

-0.51* 

0.01 

-0.72* 

-0.43* 

2 ^stai- 

-0.53* 

0.09 

-0.41* 

0.12 

-0.37* 

0.22 

-0.58* 

-0.20 


positions whose 5F£’inoi(L12) is well below the mean. They have an abnormally high cirrus correctiorp^in Equation A1 


and their Ssfr(L12) is also the lowest in our sample, < 2 x 10 ^ Mq Myr ' pc These SFR leve ls are similar to the 


< 1 X 10 ^ Mq Myr ‘ pc ^ threshold below which hybrid SFR estimators may become unreliable (L12i. On the other 


hand, there are a few positions whose 5T'£’nioi(L12) is well above the mean. They have some of the highest HQ'-to-24 jim 
intensity ratios in our sample. We discuss them in the next point. 

• The Ha-bright positions mentioned above tilt the otherwise flat fit to S FFmoKHa). Our sample does not include enough 
Ha-bright positions to assess a possible difference in S FE^oi between the two SFR regimes dominated by obscured and 
unobscured star formation. In particular, the unobscured star formation only amounts to ~ 5% - 55% (median 30%) of the 
total SFR(L12) in positions detected in HCN. In our case, pure IR and the L12 values likely provide more reliable estimates 
of the total SFR. Furthermore, in their thorough study of HERACLES galaxies |L 13 1 found that the environmental trends in 
the molecular depletion time (i.e., the inverse of SFE^oi) were not significantly sensitive to the chosen SFR tracer. 

Something similar is found when payi ng attention to the rank coefficients between the star formation efficiencies and the ISM 
parameters, which are listed in TabIe|A2| On the one hand, S FEdense remains strongly and significantly correlated with Smoi/Satom 
and Sstar in all cases On the other hand, S FE^oi remains uncorrelated with the ISM parameters, except when Ha is chosen as 
SFR tracer. 
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Their cirrus emissivity is normal, but the Sdiist/^ 24 /im ratio is high. 
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B. DATA 


We summarize here the results of our observations. Figs. |Blf[B4| show the location of the observed posi tions on th e HERACLES 
CO(2-l) maps of the target galaxies. The HCN(l-O) and CO(l-O) final spectra are shown in Eigs. B^B19 The velocity- 
integrated intensities of the two lines are listed in Table |Bl[ 
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Table B1 

HCN(l-O) and CO(l-O) velocity-integrated intensities 


Name 

NGC 

off 

r25 

2hcn 

(K km s“') 

2 coio 

(Kkm s“*) 

0628 

1 

0.00 

0.17 

(0.03) 

7.5 

(0.32) 


2 

0.27 

0.09 

(0.03) 

5.8 

(0.26) 

2146 

1 

0.02 

4.44 

(0.13) 

101.2 

(1.45) 


2 

0.42 

<0.24 


12.1 

(0.95) 

2403 

1 

0.06 

<0.12 


7.5 

(1.58) 

2798 

1 

0.20 

1.23 

(0.11) 

31.3 

(1.54) 

2903 

1 

0.01 

3.13 

(0.10) 

54.1 

(2.01) 


2 

0.14 

0.77 

(0.06) 

40.0 

(1.67) 


3 

0.14 

0.52 

(0.04) 

32.0 

(1.79) 


4 

0.37 

<0.12 


5.5 

(1.11) 

2976 

1 

0.39 

0.15 

(0.03) 

9.6 

(0.22) 


2 

0.40 

0.10 

(0.03) 

7.3 

(0.24) 

3034 

1 

0.02 

29.55 

(0.18) 

507.1 

(1.51) 


2 

1.76 

<0.16 


9.0 

(0.98) 

3049 

1 

0.04 

0.11 

(0.03) 

<4.3 


3077 

1 

0.01 

<0.12 


4.9 

(0.78) 

3184 

1 

0.32 

<0.08 


6.3 

(0.12) 


2 

0.25 

<0.16 


5.2 

(0.53) 


3 

0.01 

0.19 

(0.03) 

9.2 

(0.12) 

3198 

1 

0.01 

0.19 

(0.04) 

12.3 

(0.69) 

3351 

1 

0.01 

2.43 

(0.09) 

42.8 

(1.13) 

3521 

1 

0.12 

0.59 

(0.08) 

35.2 

(0.81) 


2 

0.11 

0.76 

(0.11) 

39.0 

(0.79) 


3 

0.41 

0.15 

(0.03) 

12.3 

(0.41) 


4 

0.38 

0.28 

(0.05) 

13.7 

(0.42) 

3627 

1 

0.01 

2.42 

(0.11) 

58.2 

(0.78) 


2 

0.19 

1.59 

(0.08) 

53.8 

(0.64) 


3 

0.51 

0.30 

(0.04) 

14.3 

(0.46) 

3938 

1 

0.04 

0.28 

(0.05) 

7.6 

(0.45) 

4254 

1 

0.44 

0.21 

(0.03) 

9.4 

(0.28) 


2 

0.16 

0.77 

(0.06) 

24.0 

(0.40) 


3 

0.86 

<0.12 


4.0 

(0.20) 

4321 

1 

0.45 

0.11 

(0.03) 

9.6 

(0.33) 


2 

0.59 

0.11 

(0.03) 

5.2 

(0.23) 


4 

0.01 

2.95 

(0.11) 

55.9 

(0.73) 


5 

0.29 

0.41 

(0.04) 

12.3 

(0.36) 

4536 

1 

0.02 

1.45 

(0.10) 

57.4 

(0.85) 

4569 

1 

0.01 

2.57 

(0.13) 

73.9 

(0.79) 

4579 

1 

0.01 

1.20 

(0.07) 

18.6 

(0.74) 


2 

0.48 

<0.16 


6.9 

(0.49) 

4631 

1 

0.24 

0.84 

(0.10) 

47.3 

(0.64) 


2 

0.67 

<0.12 


8.1 

(0.39) 

4725 

1 

0.06 

<0.24 


6.1 

(0.62) 

4736 

1 

0.00 

1.16 

(0.13) 

35.1 

(0.74) 


2 

0.18 

0.31 

(0.09) 

16.2 

(0.45) 

5055 

1 

0.00 

2.21 

(0.11) 

67.4 

(0.72) 


2 

0.18 

0.33 

(0.06) 

16.9 

(0.44) 


3 

0.32 

0.22 

(0.04) 

15.7 

(0.34) 


4 

0.42 

<0.12 


8.5 

(0.31) 


5 

0.48 

<0.12 


12.3 

(0.36) 

5194 

1 

0.15 

2.41 

(0.09) 




2 

0.50 

0.63 

(0.06) 




3 

0.00 

4.88 

(0.10) 



5457 

1 

0.00 

0.64 

(0.03) 



5713 

1 

0.11 

1.22 

(0.06) 

39.2 

(0.64) 

6946 

1 

0.00 

9.57 

(0.22) 

224.3 

(0.87) 


2 

0.42 

0.47 

(0.05) 

23.5 

(0.22) 


3 

0.73 

<0.15 


2.5 

(0.20) 


4 

0.47 

0.58 

(0.05) 

23.7 

(0.24) 


5 

0.75 

0.17 

(0.02) 

9.6 

(0.23) 

7331 

1 

0.15 

1.16 

(0.11) 

40.5 

(0.35) 


2 

0.13 

0.93 

(0.10) 

39.9 

(0.50) 


Note. — Errors at Icr-level are written within brackets. 
We give 4cr values for non-detections. We omit position 3 
in NGC 4321, which was only observed in CO(l-O). 
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Figure B2. Same as Fig.j^forNGC 3034, NGC 3049, NGC 3077, NGC 3184, NGC 3198, and NGC 3351. 
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Figure B3. Same as Fig.j^for NGC 3521, NGC 3627, NGC 3938, NGC 4254, NGC 4321, and NGC 4536. 
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Figure B4. Same as Fig.j^for NGC 4569, NGC 4579, NGC 4631, NGC 4725, and NGC 4736. 
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Figure BS. Same as Fig.j^for NGC 5055, NGC 5194, NGC 5457, NGC 5713, NGC 6946, and NGC 7331. 
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Figure B6. HCN(l-O) (left) and CO(l-O) (right) spectra at the positions observed in NGC 0628, NGC 2146, and NGC2403. The vertical dashed lines delimit 
the line windows used to reduce the data and calculate the line parameters. In each panel, the red spectrum is the HERACLES CO(2-l) line at the same position, 
convolved to matched spectral and spatial resolution. The CO(2-l) line is rescaled so that its velocity-integrated intensity within the line window is: (1) the same 
as that of the other line, if the latter has SNR\ine ^ 4; (2) equivalent to a 4cr value for the other line, otherwise. 
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Figure B7. Same as Fig.j^for NGC 2798 and NGC 2903. 
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Figure B8. Same as Fig. [B6lfor NGC 2976, NGC 3034, and NGC3049. 
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Figure B9. Same as Fig. [B6lfor NGC 3077, NGC 3184, and NGC 3198. 
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Figure BIO. Same as Fig.[^for NGC 3351 and NGC3521. 
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Figure Bll. Same as Fig.[^for NGC 3627 and NGC 3938. 
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Figure B12. Same as Fig. |B6| for NGC 4254. 
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Figure B13. Same as Fig. [B6lfor NGC 4321 and NGC 4536. 
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Figure B14. Same as Fig.[^for NGC 4569, NGC 4579, and NGC 4631. 
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Figure B15. Same as Fig. [B6lfor NGC 4725 and NGC 4736. 
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Figure B16. Same as Fig. [B6lfor NGC 5055. 
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Figure B17. Same as Fig. [B6lfor NGC 5194 and NGC 5457, and NGC 5713. 
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Figure B18. Same as Fig. |B6| for NGC 6946. 
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Figure B19. Same as Fig. [B6lforNGC 7331. 










































































































